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Abstract. Within the framework of field-theoretical description of second-order 
phase transitions via the 3-dimensional O(N) vector model, accurate predictions for 
critical exponents can be obtained from (rcsummation of) the perturbative series 
of Renormalization-Group functions, which are in turn derived — following Parisi's 
approach — from the expansions of appropriate field correlators evaluated at zero 
external momenta. 

Such a technique was fully exploited 30 years ago in two seminal works of Baker, 
Nickel, Green and Meiron (PP"E1)> which lead to the knowledge of the /3-function up 
to the 6-loop level; they succeeded in obtaining a precise numerical evaluation of all 
needed Feynman amplitudes in momentum space by lowering the dimensionalities of 
each integration with a cleverly arranged set of computational simplifications. In 
fact, extending this computation is not straightforward, due both to the factorial 
proliferation of relevant diagrams and the increasing dimensionality of their associated 
integrals; in any case, this task can be reasonably carried on only in the framework of 
an automated environment. 

On the road towards the creation of such an environment, we here show how a 
strategy closely inspired by that of Nickel and coworkers can be stated in algorithmic 
form, and successfully implemented on the computer. As an application, we plot the 
minimized distributions of residual integrations for the sets of diagrams needed to 
obtain RG-functions to the full 7-loop level; they represent a good evaluation of the 
computational effort which will be required to improve the currently available estimates 
of critical exponents. 
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1. Motivation 

A great achievement of the Renormalization-Group (RG) approach to critical 
phenomena (introduced in [3], jlj, [5], [Zj, [H]) is the realization that the universal 
behavior of many different physical systems can be explained by the existence of a 
common large-distance fixed point — or infrared, or IR fixed point — of RG-equations; 
these equations describe the evolution of the effective Hamiltonians for such systems 
when the short- distance degrees of freedom — also called ultraviolet, or UV — are 
increasingly summed up. In particular, it turns out from this idea that a quantitative 
description of physical universal quantites such as critical exponents and critical 
amplitudes can be obtained by considering the RG-evolution of a euclidean quantum 
field theory appropriately chosen to stay in the same universality class of the systems 
to be described. In this perspective, second-order phase transitions of systems with 
N phenomenological scalar components in the sense of Landau theory can be obtained 
by studying the IR fixed point (s) of the RG-equations for quantum field theories of iV 
interacting scalar fields. In the case of 0(iV)-symmetric interactions of (0 2 ) type, the 
IR fixed point is known as Wilson-Fisher fixed point (see p]). 

Following this line of thought, during the Cargese summer school in 1973, jHI], Parisi 
made the seminal observation that a quantitative description of the Wilson-Fisher fixed 
point can be obtained by analyzing the Callan-Symanzik RG-equations for a system of 
3-dimensional massive scalar bosons with an 0(A r )-invariant interaction of (0 2 ) 2 type, 
the so-called O(N) vector model (first published notes of this idea can be found in [TT|). 
Such a theory is a super-renormalizable one, i.e. the number of one-particle irreducible 
(OPI) diagrams with new primitive divergences is finite, and renormalization process is 
easier if compared to that required by other renormalizable theories. 

A price to pay for the simplicity of Parisi's approach is the fact that the value of 
the coupling is nonperturbatively large at the fixed point; thus, the perturbative series 
obtained in this framework have to be resummed by means of a suitable method (the 
one usually chosen being the so-called Borel summation) . We do not dive here into the 
details of these resummation techniques; we only mention that for this field-theoretical 
model it has been rigorously proved in that resummed perturbative series converge 
when their order increases; thus, the precision of the estimates obtained within Parisi's 
approach is only limited by our practical ability in generating and parametrizing the 
necessary amplitudes up to the desired order, and by the precision we can reach when 
computing the high-dimensional integrals in terms of which amplitudes are expressed. 
As an additional remark, one should bear in mind that, to avoid instabilities in the 
resummation process at a given perturbative order, the numerical precision of each 
term of the series should be greater than that desired for the final result, and greater 
than that of subsequent terms as well4 

J Being the subject of critical phenomena so rich and complex, it is clearly impossible for us to give 
in this context much more than a very limited introduction to it; however, the interested reader can 
easily find complete and pedagogical expositions in a variety of textbooks (see e.g. Eli E3> 
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The panorama of results in the computation of RG-functions for the 3- dimensional 
O(N) vector model is dominated by the work of B. G. Nickel and collaborators. In 1976- 
1978 — see pQ and [2] — , making use of a large and cleverly- arranged set of remarkable 
computational simplifications partially described in [20], Baker, Nickel, Green and 
Meiron were able to push the calculation from the 2- loop- (see jTHJ and [21]) up to 
the 6-loop- level. (See also [22] for a compilation of diagrams, weights and associated 
amplitude values; [23] is an application of the same type of data to a different context.) 
In 1991 Murray and Nickel [21] completed the evaluation of an improved estimate of 
the anomalous dimensions rj and r] 2 including partial 7-loop results, but the full 7- 
loop computation of RG /3-function and critical exponents is still an open problem. It 
should be noticed that from a technical point of view — and in particular when taking 
into account the lack of available computer power at the time — the 6-loop calculation is 
already an outstanding masterpiece in its own right: in fact, it requires the high-precision 
numerical evaluation (with 8-10 significant digits) of a thousand diagrams, which were 
apparently parametrized by hand as integrals with up to 6 residual dimensions, in the 
careful attempt of finding the most advantageous representation for each amplitude. 

As a matter of fact, a very conspicuous body of literature relies on [2], [21] and the 
compilation [22] for the calculation of critical exponents as well as of other quantities 
(see for example j2S] and [2H], [2*7] . [2*5]). However, at least in our knowledge, no other 
group has been able to reproduce these computations independently so far: this fact 
has been our main motivation for resuming the problem, in the hope of being able to 
verify — and possibly extend — the results available at present. 

The high number of diagrams which must be evaluated to complete the calculation 
of the 7th loop (~ 4000, see section I3.2.2|) makes it immediately clear that the only 
realistic hope to push forward our knowledge relies on finding a way of automating 
computations by means of computer techniques; the goal of this paper is then to 
perform the first steps towards the build-up of such an automated environment, both 
for generation and parametrization of all needed Feynman amplitudes. Article structure 
is as follows: 

• In section [2] we introduce all relevant information about the 3-dimensional O(N) 
vector model, together with the renormalization scheme and the RG-equations our 
computations are based on 

• In section E] we briefly discuss how to generate all needed Feynman diagrams, 
together with their combinatorial and O(N) factors; after that, we analyze in great 
detail the set of analytic tricks which make possible an efficient parametrization of 
the associated amplitudes. A number of practical examples is given 

• In section 0] we provide a precise algorithmic formulation for each of the tasks 
required to actually find the set of optimal parametrizations at some given 

[T7j : recent reviews are presented in ^S] and EH])- ln these texts, many other subjects relevant to 
this article — like euclidean quantum field theories, theory of Renormalization-Group and resummation 
techniques — are also extensively covered. 
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perturbative level 

• In section E] we describe our computer implementation of such techniques, and the 
results we obtained so far. A critical discussion follows about how the set of chosen 
simplifications influences both the complexity of the parametrization code and the 
minimality of the resulting set of parametrizations. 



2. The 3-dimensional O(N) vector model 

2.1. Generalities 

The 3-dimensional O(N) vector model is described by a classical (bare) action which is 
invariant under a transformation of the scalar fields 0' B (i = 1, . . . , N) according to a 
fundamental representation of O(N): 

Wb) := / d*x Q(V0 B ) 2 + ^ + ^) 2 ) . (1) 

The quantization of the model is obtained via functional integration, by considering the 
(bare) generating functional, 

Z B (J B ) := J [d4> B ]exp (-S B (<f) B ) + J J B <p}j . (2) 

An appropriate regularization — parametrized by some ultraviolet cutoff denoted with 
A — is necessary to give meaning to functional integration. We assume that such a 
regularization exists, but to maintain generality and simplicity of notations we will 
not show its details here. As well known, perturbation theory for the bare model is 
plagued by ultraviolet divergences in the limit A — > oo, and must be complemented by 
renormalization; that amounts to a reparametrization of bare (ultraviolet) quantities in 
terms of renormalized (infrared) ones. In this section we content ourselves with briefly 
recalling only some essential points about the renormalization of the O(N) vector model; 
all other necessary definitions and some additional technical details about it, including 
a derivation of Callan-Symanzik equations in a generic massive renormalization scheme, 
are reported for reference in | Appendix A| 

Applying standard power counting techniques to Feynman diagrams obtained from 
the perturbative expansion of the O(N) vector model, it is easy to compute the 
superficial degree of divergence to for the Feynman amplitude A® associated to a 
connected one-particle irreducible, or OPI, diagram D with E external legs, V 2 insertions 
of operator J 2 and V4 vertices with four legs. The obtained expression reads 

u = 3 - | - V 4 - 2V 2 ; (3) 

roughly speaking, u parametrizes the behaviour of the amplitude in the limit A —>■ 00, 
where A® ~ A w (see for a more rigorous definition). It turns out from (jSJ) that the 
model is super-renormalizable, i.e. it has only a finite number of superficially divergent 
amplitudes. When restricting to the conditions E > — as it is always the case in this 
article — and L > 0, amplitudes must satisfy E — 2, V2 — and V4 = 1, 2 to have u > 
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and to possibly be superficially divergent; these constraints are verified in our model 
only by the three graphs listed below: 




which we will call in the rest of this article — from left to right — "tadpole", "cactus", 
and "sunset", respectively. It should also be noticed that these divergent graphs require 
a mass renormalization, but no divergent contribution affects the derivative of the two- 
point function w.r.t. p 2 , i.e. wave-function renormalization is not needed to make the 
theory finite. 



2.2. Intermediate scheme 

From a practical point of view, the most effective way of proceeding is to first choose 
an intermediate renormalization scheme, with the goal of making the theory finite and 
the calculations as simple as possible; only after that — in the end — one will switch 
to the final renormalization scheme presented in section 12.31 More in detail, our own 
intermediate scheme (labeled with suffix i and referred to as "I-scheme") is defined by 
the following conditions: 

0B = 0! (4) 

B =[0 2 L (5) 

A B = m^j (6) 

m\ = m\ + 5m 2 ; (7) 
the mass counterterm is in turn defined as: 

K^-(-Q- + -e-|J- (8) 

From such definition it follows that in the I-scheme the renormalized contribution of the 
"tadpole" fully vanishes, while the renormalized contribution of the "sunset" vanishes 
at p = 0. As a consequence, the renormalized contribution of the "cactus" diagram 
vanishes too, because the "cactus" amplitude factors in terms of a "tadpole" amplitude 
times a finite term. (Please remark that in (jHJ) we did not write the explicit form of 
5m 2 , being it dependent on the choice of regularization.) 

As a first important property of our intermediate scheme, an inspection of 
counterterms listed in (jHJ) readily shows the absence of potential overlapping divergent 
(sub)graphs. Thus, the procedure of renormalization amounts in the I-scheme to 
mechanically replacing all the instances of divergent graphs with their renormalized 
counterparts when they occur as subgraphs embedded in larger graphs; in practice, all 
needed operations can be performed by just putting "tadpole" subdiagrams to zero, and 
by replacing all occurrencies of the "sunset" with a renormalized version, obtained from 
the original one by subtracting the value of the diagram in zero. 

This trivialization of renormalization is a fondamental assumption in the automated 
framework presented in this paper, and its consistent practical advantages will be silently 
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used everywhere throughout all our parametrization algorithms; it should be noticed 
that this property — which is far from being always guaranteed — depends in general 
both on the structure of primitive divergences of the model considered and on the type 
of subtractions that are performed in a given renormalization scheme§. 

Another remarkable property of the I-scheme is that <\> x and [<f) 2 ]i do not renormalize, 
i.e. Zfo = Z^ 2l = 1; this fact implies the vanishing of corresponding anomalous 
dimensions 

Vi(9i) = VM = (9) 

(see |Appendix A| for general definitions of Z-terms and anomalous dimensions). 
Triviality of Z-terms implies in turn the following set of relations, valid for E > 0: 



T[ Ey2) (p; m 2 , g,) = T^'^ (p; qr; m 2 B = m\ + 5m 2 , A B = m^) . 



,(-E,V 2 ) 



(10) 



Remark in the equation above the use of p, q to denote respectively the collections 
(p 1; . . . , p E ), (qr 1 , ... , qy 2 ); please consult Appendix A| for more details on our 
notations. 

The Callan-Symanzik operator in the I-scheme, T> z , is defined by specializing to 
such a scheme the general definition in (|A.19J) : 

d 



dm l 



A B ,A 



Applying T* l to both sides of (JHJ) the expression for the /^-function is easily derived: 

(3 1 (g 1 ) = V 1 [g 1 ] = -g 1 . (12) 
The exact Callan-Symanzik equation in the I-scheme (assuming E > 0) then reads: 



d 







m, 



dm l ^ dg x 



if™ (p; q) = mfoCfc) rf >* +1 > (p; q, q Vi+1 = 0) (13) 



(E,V 2 +1) , 



where the cx-function is obtained by applying its definition (jA.26|) to (JJJ), 



(<?i) = 2 + 



"Pi 



AT + 2 2 
^i- 



(14) 



24vr " 288vr 5 

(It should be noticed that, regardeless of the dependence of 5m 2 on the cutoff A, the 
expression of o Y (g^ does not depend on the choice of regularization in the limit A — > oo, 
as it should be for all RG-functions in a well-behaved renormalization scheme.) 

In spite of their apparent simplicity, Callan-Symanzik equations in the I-scheme are 
not at all trivial! For instance, specializing to (E, V2) = (2, 0) at p = and introducing 
the rescaled function 

«^-^(0;m Ilft )K 



we obtain the exact relation 

d 



gi 



t[ 2 > \ 9i ) = aM ir"(0;« = 0;a) 



(2,1), 



(15) 



(16) 



§ See e.g. |30] for a pedagogical discussion about overlapping divergences and Zimmermann's forests 
formula. 
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which can be used to link two quantities involved in the computation of RG-functions, 
see discussion below. 

Before closing this section some comments are in order to motivate our own choice 
of the I-scheme as intermediate renormalization scheme. 

It should be clear from the considerations presented above — and in particular from 
the definition of mass counterterm in (jSJ) — that the main peculiarity of the I-scheme 
is that a minimal number of counterterms is introduced. This property results in a 
minimal amount of renormalized diagrams and turns out to be invaluable to simplify the 
automatization of the parametrization process, expecially in view of cost optimization 
- see section 0] This simplicity, together with the desire of testing existing results 
in the most independent way as possible, mainly motivates our choice of the I-scheme, 
which is in fact quite different from the one employed in j22] . 

The price we have to pay for this conceptual simplification is that in this scheme 
we may have more diagrams to evaluate than in other possible schemes; in our case, 
for example, the renormalized high-order contributions to r[ 2,0 ' ) (p = 0) are not trivially 
vanishing — apart from the "sunset" graph and all graphs with "tadpole" insertion(s). 
In particular, one relevant difference w.r.t. the scheme used in [22] is that in our I- 
scheme 2-point connected OPI (sub-)graphs with external lines connecting to the same 
vertex (which we refer to as to generalized tadpoles ||) give in general rise to non-zero 
contributions, and must thus be evaluated. However, this complication is not a dramatic 
problem for at least two reasons: 

(i) the expression of r[ 2 '°\p = 0) can actually be obtained from that of If' ^(p = q = 
0) by use of the exact Callan-Symanzik equation in I-scheme, (fTo^) 

(ii) anticipating section 12 .3\ we remark that in addition to that of r| 2,0 ' ) (p = 0) also 
the expressions for r^ 4 , rp' 1 '' and dv[ 2 ' ^ /dp 2 — all evaluated at zero external 
momenta — are required to obtain the desired RG-functions in the massive scheme. 
It turns out in practice that the actual complexity of the evaluation of r| 2 ' ^ (p = 0) 
is relatively small if compared to the computational costs of the other three needed 
correlators just mentioned. 

Motivated by those considerations, we will include rj 2,0 ' ) in our cost analysis of section 
planning to use (fTT)|) as a consistency check on future results. 

2.3. Parisi's massive scheme 

The standard massive scheme — suffix m — , as introduced in [IT], is defined by the 
following normalization conditions: 



p=0 

|| They are called "Hartree-type self-energy insertions" in |22j . 




(17) 
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rL ' ^(p'i m Mt 9mj 



p=0 



-n(2,l)/ 2 a 

Tm (p, -p; q;m M ,g M/ 



m M g M (19) 
= 1. (20) 

p=q=0 

(Please remark that we are adopting here the same conventions introduced for bare 
correlators in equations (jA.7HA.9p .) The relation with bare quantities, in the case E > 0, 
follows from (jA.18j) : 

v 2 



J <f),M) 



,B/2 



^12 



(21) 



Callan-Symanzik equations (jA.21|) . and definitions of RG-functions as in equations 
(|A.23HA.26|) . hold in this scheme with the suffix replacement r — > m. Specializing ()A.21j) 
to (|17|) we get an additional relation specific to this renormalization scheme: 



2 - 77m(#m) = Cm(S'm)- 
The other RG-functions are obtained as follows: 
(i) using (fTU|) . ()21|) and normalization conditions 



ar (2,0) 

^-Cp, -p. 



(22) 



we obtain the relations 



r{ 2,1) (oj o) 



p=0, 



</> 2 ,M 
fl'M = 



r (4,0) (Q) 



(23) 

(24) 
(25) 



r (2,0) (Q) 



(2,0) 



m: 



dp 

r (2,0) (Q) 



— IP. "Pi 



p=0 



ar (2,0) 
"'" 1 ^2-(P'-P) 



(26) 



p=0 



(ii) by definition (and use of chain rule) we observe that 

9g M 



dgi 



Vm^m) = V M [\ogZ, 



<P,m\ 



d log Zfa 
dgi 



V u [log Zf 



dlogZ* AA 



(27) 
(28) 

(29) 



(iii) the expression for V M [g^ is readily derived by applying V M to both sides of relation 



A B = rriigi, and using the identity 
V-JmA . 1 



log- 



1 dlog- 

1 " 2 ^ -oi 



(30) 
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As a result we obtain: 




-i 



Z> M [ft] 



9i\ 1 



(31) 



Equation (}3*Tj) together with equations (j23H26j) and equations (I27H29|) allow to determine 
RG-functions in the massive scheme from the knowledge of r[ 2 ' \ of its derivative w.r.t. 
p 2 , of Ti 4,0 '* and rj 2,1 '*, all evaluated at zero external momenta. 

To compute critical exponents one must first resum with some appropriate method 
(e.g., Borel summation) the divergent perturbative series obtained for (3 M (g u ), for 
anomalous dimensions r) M (g u ), t^m^m), and — as a typical check — for other series 
derived from anomalous dimensions by usual scaling and hyper-scaling relations among 
exponents, like e.g. v m {.9m) = 1/(2 + ri2 M (g M ) — r) M (g M )). Critical exponents are 
then obtained by evaluating the resummed series at g M = g^, the non trivial zero 
of Pm{9m), which corresponds to the critical region. In spite of the fact that they can 
be obtained from scheme-dependent quantities, exact critical exponents are universal 
numbers. (Notice nevertheless that the speed of convergence of approximations obtained 
from resummation of perturbative series at fixed finite order might vary by choosing 
different renormalization schemes.) 

3. Setting up the problem 

3.1. Flowchart 

As shown in previous section, performing the calculation of RG-functions is tantamount 
to computing perturbative series up to the desired order for correlators r| 2 '°' ) , rj 4,0 -*, 
If and <9r[ 2 '°' ) /dp 2 , all evaluated at zero external momenta. 

In figure ITTI we present a flowchart where all the tasks contributing to the structure 
of our automatic framework are mentioned. Quite schematically, we can identify three 
main steps: the generation of all needed diagrams with their associated weights and 
O(N) tensor structure, the parametrization of amplitudes associated to diagrams and, 
finally, the numerical integration of generated code. 

The problem of generating diagrams and combinatorial factors is computationally 
hard but standard, and will thus be addressed only briefly at the beginning of section 
El the last part of section 01 is devoted to prepare the theoretical ground for section El 
where we present the algorithmic implementation of an efficient strategy to find for each 
amplitude integral representations with low associated dimensionality. As anticipated, 
in this paper we will not deal at all with the challenging problem posed by the numerical 
integration of obtained parametrizations, leaving this task to future publications. 

3.2. Graphs, weights and tensors 



3.2.1. Representation of graphs While many good computer packages for manipulating 
graphs do exist (see for instance (SI), some of them being specifically tailored to 
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Other factors 
of f actorizable 




First factor of factorizable 
i Not factorizable 



Combinatorial weights 




O(N) tensors 





dp 2 



Filter for chain diagrams 



Factorizable 



Lookup at lower order 




Not factorizable 



Generation-related Tasks 



Renormalization 



List of known vertex functions 



2 



Substitutions 



Parametrization-related Tasks 



Choice of momenta 

I 



Minimization of cost 



Choice of angular measure 

z 



Emission of C++ code 



i 

Numerical integration 



Figure 1. Flowchart of our automatic framework. Its details are clarified in sections 

nrn 



applications in theoretical physics (one such example is [32]). we have anyway decided to 
re-implement our own graph tool from scratch. The main motivation for such an effort 
can be found in the fact that, since our entire parametrization chain is formulated in 
terms of graphs, we needed much more extended capabilities than those usually found 
in usual Feynman-graph manipulation tools; furthermore, we badly needed efficiency, 
because some of the stages of parametrization production are computationally very 
expensive (see section IB. 6. II and the end of section fOj) . 

All these considerations brought us to the conclusion that only a new library, written 
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in terms of an efficient low-level programming language and providing a set of graph- 
theoretical operations specifically tailored to our needs, could solve our problems — in 
the end, C++ has been our final choice of language. 

Here are the basic design choices which have been formulated for our library: 

• it allows graphs with coloured vertices and coloured lines to be described (in our 
situation, coloured vertices are used to model external lines and effective vertices 
— see section l4~21 — , while coloured lines are needed i.e. for symmetries of arguments 
of effective vertices and cosine diagrams, see section 14 .4|) 

• it is based on an adjacency-matrix representation (we briefly recall that for diagrams 
where only a single kind of lines is allowed, each entry of the adjacency matrix 
A v y gives the number of lines joining vertex v and vertex v'; in more complicated 
situations A v y encodes both the number and the colours of such lines) 

• it knows the concept of canonical representative for a graph. The basic idea here 
is the fact that the same graph can have in general different representations, 
the one obtained from the other by a relabeling of graph vertices; since such 
a relabeling defines an equivalence relation, the canonical representative is - 
as usual in similar situations — the special representative picked up to label all 
the equivalent^ elements belonging to the same class. It should be noticed that 
implementing the concept of canonical representative was particularly important in 
our context: in fact, the definition of the representative can be used to provide an 
ordering relation for graphs, which is in turn essential to a good cooperation of our 
library with the data-types and the facilities offered by C++ standard library |33j . 
However, while the knowledge of the canonical representative for a graph gives 
many advantages, it has the main drawback that its computation can turn out to 
be very expensive (in fact, this operation can imply the inspection of a number 
of graphs up to the factorial of the number of vertices, such being the maximal 
number of possible vertex relabelings) . 

On the other hand, many algorithms exist which allow to reduce the cost of finding 
the canonical representative (even if each algorithm known so far suffers from some 
"difficult" graphs); these algorithms usually consist in "colouring" vertices on the 
basis of some vertex property which is invariant by vertex relabeling, so to partition 
the vertex set into smaller subsets, and to reduce — sometimes drastically — the 
number of permutations to be checked. Our personal choice among the many 
algorithms proposed in the literature (see for example jMj, jUSj and [HE] for an 
overview) has been to implement a twofold colouring, based both on leg partitions 
and on minimal distances (refer to [34J for a description of the latter); these criteria, 
though being not as efficient for large graphs as those used in [3*5] . are far less 
complicated from the point of view of software implementation, and nonetheless 
quite efficient for the small Feynman diagrams we need to cope with. 

Some of the graph-theoretical operations provided by our library are for example the 

If Or, following standard graph terminology, isomorphic. 



Automatic parametrization of amplitudes 



12 



Table 1. Number of diagrams (OPI and without tadpoles) contributing to the two- 
point and four-point function up to 8 loops in perturbative order. 











Incremental number 




Total number 


Loop number 





1 


2 


3 4 5 6 7 


8 


6 7 8 


p(2,0) 








1 


2 6 19 75 317 


1622 


103 420 2042 


r (4,0) 


1 


1 


2 


8 27 129 660 3986 


26540 


828 4814 31354 



identification of connected components, the computation of the order of automorphism 
group, the computation of spanning trees, the enumeration of all occurrences of a given 
subgraph in a larger graph (presented as an ordered sequence, that is as an iterator in 
C++ terminology, see [HB]), and the enumeration of all cycles, of all cycle bases and of all 
paths joining two vertices (in the form of iterators as well). All these capabilities will 
be essential to implement algorithms presented in section 0J 

3.2.2. Generation of graphs The task of generating all inequivalent graphs which 
satisfy a predefined set of topological conditions is known since a long time to be a 
very hard one, due to the fact that all known generation algorithms are also plagued 
by the problem of isomorphic copies: diagrams corresponding to the same canonical 
representative are in general produced more than once, leading to a worser and worser 
inefficiency of the process as the number of vertices of the graph increases. Various 
methods have been proposed in the literature; here again, some packages are available 
(see for example [SZ|)- In the spirit of jHH], we generate the set of diagrams in a non- 
recursive way, taking advantage of our knowledge of the canonical representative to 
prune a large number of isomorphic copies out of the generated tree (see [HI]). The 
performance of such algorithm is very satisfactory, at least for the typical number of 
vertices we are interested in. 

In table [U we quote some results about the number of diagrams needed to compute 
the two-point and four-point functions. The main facts we can deduce from these 
numbers are the following: 

(i) to complete the evaluation of the RG /5-function at 7-loop level 3986 more diagrams 
must be evaluated just for the four-point function vf^ alone 

(ii) improving the computation of anomalous dimensions r] M and r] 2M in terms of g l to 
8 loops would require the evaluation in p = of both values and derivatives w.r.t. 
p 2 of the 1622 diagrams contributing to r| 2,0 ' ) , plus — eventually — the evaluation 
of diagrams contributing to T\ \ (p = q = 0). (For a more precise statement 
of the problem, please refer to the discussion at the end of section 12.21 ) These 
diagrams have more lines and, accordingly, are more difficult to evaluate than those 
contributing to the seven-loop four-point function, see section ED 

In any case, results in table Estate clearly that neither of the two tasks presented above 
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can be performed without the help of an automated framework, able to supply the 
parametrization and the evaluation of all needed amplitudes. 

3.2.3. Symmetry factors When generating Feynman diagrams for the perturbative 
expansions of correlators in quantum field theories, one is interested in describing graphs 
whith unlabeled internal vertices, and unlabeled lines as well, because graphs differing 
by such relabelings contribute with the same amplitude; this causes the problem of 
computing the symmetry factor for a given Feynman graph, which accounts for the 
multiplicity of such identical contributions. 

Practical computation of symmetry factors can be very difficult, since it requires 
the knowledge of the order of the automorphism group of the graph of interest; in our 
case, symmetry factors are derived directly from the adjacency-matrix of the graph, 
following the directions given in |4Uj . 

3.2.4- O (N)- factors The computation of traces for O(N) group is quite simple, and 
can be carried out in many ways; for example, tensors can be represented as graphs, 
and their contractions can be performed in full graphical form as well. Since this point 
does not pose particular problems, we will not insist on it here. 

3.2.5. Consistency checks Of course, all possible precautions must be taken against 
the unpleasant possibility of accidentally omitting a graph, or computing a wrong 
symmetry/O (N) factor. The consistency of our results has been carefully checked in 
many ways; among them we recall the fact that from zero- dimensional field theory some 
sum rules can be deduced, allowing to test both combinatorial and 0(N) factors (see 
[Ti] for a description of this technique). As an aside, we notice that the introduction 
of a mass counterterm in such a context allowed us to directly obtain sum rules for 
diagrams without tadpoles (and an analogous solution has apparently been adopted in 

E21). 

3.3. Which parametrization to choose? 

Ideally, in our approach to the computation of RG-functions the "optimal" 
parametrization for a given Feynman amplitude is the one minimizing the CPU-time 
which is needed for its numerical evaluation, up to the desired precision. In practice, 
many factors come into play while trying to quantify in a precise way the various 
contributions to such a computational cost; for example: 

• for a given amplitude, any choice of parametrization induces its own integral 
representation, with its associated specific dimensionality 

• an estimate of the number of evaluations of a given integrand needed to compute 
the value of a multidimensional integral within a specified numerical accuracy is in 
general not precisely known a priori, since it depends in an unaccessible way on 
the regularity properties of the integrand function itself 
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• parametrizations of the same amplitude in terms of integrals of the same 
dimensionality may be obtained — see section 14.21 — by choosing in various ways 
the set of replaced effective-vertex functions; these functions may have in general 
very unequal evaluation timings 

• the same vertex function may possibly require a very different amount of floating- 
point operations when it is evaluated at different values of external momenta. 

All these considerations suggest that we should not rely on a too much refined definition 
of the cost function whose minimum characterizes the optimal parametrization; on 
the other hand, the typical asymptotic behavior of the error estimate in deterministic 
multidimensional numerical integration, 



suggests that it is the dimensionality of the integration which sets the basic scale of 
the number of required evaluations of the integrand. Thus, we will discard all other 
possible definitions, basically sticking as our cost function to the dimensionality of the 
final integral needed to compute an amplitude — even if for technical reasons the actual 
definition (given in section 14.11 and used throughout the article) will need to be slightly 
more refined, its spirit remains the same. 

Having in mind from now on the goal of minimizing the dimensionality of the 
integral representations of our amplitudes, let us evaluate such a dimensionality for four 
standard families of parametrizations: we will thus compute the basic cost of momentum 
representation [EC Momentum), that of position representation [EC position), and those of 
Schwinger (i3C schwinger) and Feynman (BCF eynma _ n ) representations. Such costs are here 
called basic to stress the fact they refer to the dimensionality of integrals as directly 
obtained from standard parametrization techniques, without considering the possible 
use of additional special simplifications (in the case of momentum representation, for 
instance, a more refined choice of angular measure could be supplied; such a choice 
will indeed play an important role in lowering the number of residual integrations, as 
described in section IP]) . 

Let us consider a connected OPI diagram D with E external lines, V internal 
vertices (V =: V4 + V2, V% and V4 being the number of vertices with 2 and 4 
legs, respectively), L loops and / internal lines. The corresponding (^-dimensional 
Feynman amplitude A® evaluated at zero external momenta and unitary masses can 
be schematically written in terms of momentum, position, Schwinger and Feynman 
representations as (respectively): 




(32) 




(33) 




(34) 




(35) 
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f i 1 

oc / Uda^iY^ai-l) (Ql(«)H, (36) 
Jo i=i i 

where we denoted by l\ the loop momenta, by ki the momentum flowing in the i-th 
internal line — which can be expressed as an appropriate linear combination of loop 
momenta — , by Pr(x) the propagator in position space, by A v y the entries of the 
adjacency matrix — i.e., the number of lines connecting vertex v to vertex v' — , and 
by Ql{&) an appropriate homogeneous function of (aj)f =1 of degree L (see e.g. |2H] for 
more details). 

Before moving to the evaluation of the basic costs for these integral representations, 
a useful property of SO (rf)-invariant integrands must be mentioned (see section E.4.21 
for a more detailed illustration in the case d — 3). When integrating a function over 
a collection of vectors v i, . . . , Vk £ the dimensionality of the original integral can 
be reduced in a standard way if the integrand function depends only on the scalar 
products among those vectors; the reduction is obtained by performing the trivial 
angular integrations which correspond to the invariance w.r.t. a global SO(d) rotation 
of integration vectors V\, . . . , vk- If K > d— 1, the number of such trivial integration is 
d(d — l)/2. Due to vanishing external momenta we evaluate amplitudes at, this property 
can be applied to momentum representation ()33j) as well as to position representation 



With this property in mind, we are now ready to complete the evaluation of basic 
costs: the dimensionality of the integral required by each representation follows from 
previous expressions, by properly taking into account the delta functions in (J3T?)) . 
and by making use of the S"0 (<i)-symmetry of (jHHJ), (jH3| (corresponding estimates (|3Tj) . 
hold for L > d — 1 and V > d, respectively) : 

^^Momentum = — (37) 

£C Position =d(V-l)-^^- = d(L+ E -*- d + V 2 ) (38) 



E - 4 

£>Cschwingcr = I = 2L -\ h V 2 (39) 

£C Fcy nman = I - 1 = 2L + + V 2 - 1] (40) 

it should be noticed that the standard topological identities 

V 4 = L + ^ (41) 
E-4 

I =2L + ^— + V 2 (42) 



have been used to eliminate V4 and / from formulas 

From expressions (|37MU|) it clearly appears that, when d > 2 and the number of 
loops L is large enough, BCp ejnman < ^Momentum ~ ^Position- In table[2]we reproduce in 
detail the costs of integrating the 3-dimensional amplitudes contributing at loop orders 
L = 6, 7, 8 to r^ 2 ' -', r 1 - 4,0 ^ and r^ 2 ' 1 -'. From this analysis of standard representations, the 
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Table 2. Basic costs (6C M omcntum, SC Posi tion, BC S 

chwingci'5 n m an 

) —for 

momentum, position, Schwinger and Feynman representations respectively — of a L- 
loop amplitude contributing to T^'^ at d = 3. 







r( 2 <°) 


r( 4 >°) 


r( 2 .!) 


L = 


6 


(15,12,11,10) 


(15,15,12,11) 


(15,15,12,11) 


L = 


7 


(18,15,13,12) 


(18,18,14,13) 


(18,18,14,13) 


L = 


8 


(21,18,15,14) 


(21,21,16,15) 


(21,21,16,15) 



apparent conclusion would be that — in the case of interest — Feynman representation 
considerably lowers the number of residual integrations if compared to momentum or 
position representations. 

However, the seminal idea — exposed in [20] — which allowed Baker, Nickel, Green 
and Meiron to obtain in pQ-[2] quite precise 6-loop estimates of the RG-functions 
consisted, in fact, in using a non-standard representation of amplitudes in momentum 
space. This new representation was obtained from the standard one by replacing in 
the amplitudes as many one-loop subintegrals as possible — corresponding to one-loop 
subdiagrams of the main diagram; in fact, since the analytic expression of one-loop 
correlators for non-exceptional momenta is analytically known from the work of Melrose 
jH] , a reduction of 3 integrations in the basic cost readily follows from each substitution 
one can make in the original amplitude. 

As a matter of fact, if we analyze in more detail the cost £>Ci mpr .Mom. of momentum 
representation when it is improved in the spirit of Nickel et al., we soon come to a 
surprising conclusion. Let us suppose that we have been able to replace in an amplitude, 
say, L rcp \ analytically known functions, each corresponding to a one-loop subdiagram. 
Assuming L — L rep i > 2, so that we can spare 3 integrations thanks to the overall 
5'0(3)-invariance of the amplitude, we then obtain the relation 

^^Impr.Mom. = 3 (L — L Tcp \ — 1) . (43) 

Now, the requirement that the obtained parametrization have no more residual integrals 
than the standard Feynman parametrization is equivalent to the following condition: 

^Clmpr.Mom. < ^Cpeynman ^ 3L rcp i > L Vj> (44) 

(notice that in the present context the lower bound on L rep i is maximized by r^ 2,0 )). More 
explicitly, when L = 7(8) the condition £>Ci mpr .Mom. < BC Feynman reads L rep i > 2(3) for 
p( 2 :0) anc [ £ rcpl > 2(2) for r^ 2,1 ) and r^ 4,0 ); this means that if we are able to perform 
enough replacements the improved momentum space representation will beat standard 
Feynman parametrization. This fact motivates the choice — Nickel's as well as ours — 
of this representation for amplitudes. 

In addition, surprises are not yet over: it turns out that in this representation - 
complemented with a convenient choice of renormalization scheme — one can naturally 
take advantage of a whole set of powerful identities and properties. For example, all the 
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occurrences of the two-loop "sunset" diagram can be replaced by the analytic expression 
of its renormalized counterpart, with a gain of 6 integrals; but much more can be 
done. The next sections are devoted to a systematic presentation of a set of additional 
simplifications and tricks which are well suited to improve this framework. In particular: 

• the idea of replacing in the amplitudes some functions known analytically — which 
we will usually refer to as to effective vertex functions, or as to effective vertices for 
short — is very powerful, and can be generalized to effective vertices other than the 
ones which were presumably used by Nickel and coworkers. A detailed analysis of 
this method and its variations is carried out in section IH.41 

• other identities valid when computing amplitudes at zero external momenta 
can be used to achieve further reductions in the final dimensionality of some 
parametrizations; they are analyzed in section 

• a good choice of angular measures — that is, of the way of writing the 
parametrization of loop integration variables in spherical coordinates — can lead 
to the analytic integration of a conspicuous number of angles. However, being 
not explicitly tied to Nickel's framework this problem is in a sense a more general 
one, and requires some additional notation which will be introduced only at the 
beginning of section HJ thus, we postpone this issue to section 14.41 

• closely related to the last point, there is the possibility of factoring and computing 
analytically some easier one- dimensional integrals of chains of propagators; this 
subject is postponed as well to section l4~4l 

However, it goes without saying that this approach to the parametrization of amplitudes 
also shows some drawbacks. 

First of all, it must be stated that we do not know exactly which set of simplification 
rules has been used in [T]-j2] to lead to such a spectacular reduction of the complexity 
of the original problem, nor we do know whether such a set can be formulated in terms 
of simple algorithmic procedures. In fact, only a hint of the basics of effective- vertex 
technique is given by Nickel and coworkers in their published literature (see [22| and 
|2()j): we have thus tried, so to say, to "reverse-engineer" Nickel's results — and preprint 
[22] m particular — to deduce or re-invent the techniques we present here, selecting 
in the end the ones which seemed to us to be the most appropriate for automation; 
however, nothing prevents some (perhaps essential) ingredients of the original approach 
from being possibly very difficult to automate, and very difficult to state in algorithmic 
form. 

Secondly, a related — and much more fundamental — problem posed by this 
otherwise appealing framework is the fact that its effectiveness cannot be proven a priori: 
no theorem exists (at least in our knowledge) stating that at a given loop order, for a 
given set of available effective vertices, some minimal number of replacements will be 
uniformly obtained for all graphs needed during the evaluation of some field-theoretical 
quantity; the risk exists that, due to stronger and stronger topological obstructions 
present in diagrams at higher orders, the set of effective vertices will prove inadequate 
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to satisfactorily reduce some "difficult" graphs. As a matter of fact, a few such graphs 
requiring more residual integrations than their relatives do appear for some choices 
of effective vertices and sets of graphs to be parametrized. In the lack of a proof a 
priori, the effectiveness of the method can be demonstrated only a posteriori with an 
explicit inspection carried out diagram by diagram; considering the very large number of 
amplitudes, the large number of different possible replacements for a given diagram and 
the tempting possibility of increasing the set of effective vertices, once again we arrive at 
the conclusion that such a complex optimization problem can be successfully analyzed 
only by means of an automatic framework, capable to handle the problem of minimizing 
the number of residual amplitudes fastly and more effectively than any human being. 
In addition, such an approach is the only one able to provide an evolutive set-up if new 
tricks are found or new parametrization strategies prove themselves to be necessary (one 
example of such a case could be an hypothetical 8-loop computation, where one could 
wish to implement, for instance, a stage choosing the best parametrization among those 
given by both momentum and Feynman representations). 

Consequently, in section 0] we will show how to build a prototype framework to 
optimize — according to a reasonable choice of simplification rules — the parametrization 
of a given amplitude, starting from a given set of user-defined effective vertices. As 
reported on in section |S] and section |Hl the knowledge of a small number of effective 
vertices already leads to very effective parametrizations for the cases L = 6 and L = 7. 

3-4- More on effective vertices 

The idea of Nickel et al. of substituting effective vertices in a given amplitude, 
thus attempting to decrease the number of residual integrations and find a cheaper 
parametrization, is very general. The underlying hypotheses are the line-locality 
of momentum representation (that is, the mapping between lines of the diagram 
and products of corresponding propagators in the integrands), integration over loop 
momenta, the triviality of renormalization (absence of overlapping divergences, implying 
that one may replace a divergent block with its renormalized counterpart) and, of course, 
the knowledge of zero-cost or — more generally — low-cost expressions for effective vertex 
functions describing subdiagrams of the original amplitude (by zero-cost expression we 
mean a subgraph which is known in terms of elementary functions and zero integrations, 
while with low-cost expression we indicate an equivalent integral representation of a 
subdiagram involving a smaller number of integrations than the one we started with). 

Various classes of low-cost vertex functions are available; they will be examined in 
the following sections. 

3.4-1- One-loop functions As already mentioned, analytic expressions in terms of 
elementary functions for one-loop functions in d = 3 are known: following our definition, 
we say that all such functions are zero-cost effective vertices. 

More in detail, in |H] one-loop correlators with E > d + 1 legs (d being the 
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Table 3. Analytic expressions for some densities at d = 3, adopting the normalizations 
of 1)45(1. In the table we used the notations M\i := rti\-\-m% and Myzi '•= mi+mz+ms; 
6(x) is the usual step function. Please remark that in the case of "sunset" diagram 
a subtracted version of ffify must be employed, while the density for the "double- 
triangle" diagram is evaluated at = 1. 
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integer spacetime dimension) are reduced to a linear combination of one-loop correlators 
with E = d legs (see also [20]). Unfortunately, due to the presence of inverse 
powers of kinematical determinants in the coefficients of the reductions, such formulas 
have a range of validity limited to non-exceptional momenta, i.e. to kinematical 
configurations such that the involved determinants are nonvanishing. One possible 
way to patch this potentially catastrophic problem is to reject such exceptional points 
during numerical integration: this strategy — suggested in [20] — has the drawback of 
requiring the compatibility of point-rejection with the chosen integration algorithm. 
As a complementary alternative, we developed in |l2] a general theoretical framework 
to deal with reductions of one-loop correlators in all kinematical situations, and we 
implemented it in an highly-optimized and robust C++ library. 

3.4-2. Functions in terms of spectral densities All higher- loop amplitudes whose 
expressions can be simplified to reduce the number of residual integrals constitute 
potentially useful effective vertices. 

A first concrete example of such a situation is obtained when considering L-loop 
amplitudes A^ 2,0 \p 2 ) contributing to T[ 2 '°\p, — p): due to their analyticity in the 
complex p 2 plane cut at p 2 < — M 2 hreshold , these amplitudes possess a dispersive one- 
dimensional integral representation, which — when renormalization is not required — 
has the standard unsubtracted form 

n^n^-l^MM) (45) 

where the density is related via Cauchy theorem to the discontinuity of the amplitude 
along the p 2 -cut, and can be obtained in terms of a sum over cuts of the diagram D by 
use of standard Cutkosky rules (see [32] for a useful introduction to such techniques). 
In table [3] we list some examples of 3-dimensional densities whose expressions can 
be obtained analytically by integration of cut-diagrams. Notice that due to its UV 
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divergence the sunset diagram must be renormalized, and satisfies in our I-scheme a 
subtracted version of ()45|) which can be obtained from it by use of the replacement 
l/{p 2 + M 2 ) -> l/{p 2 + M 2 ) - 1/M 2 . 

A second possible direction to obtain effective vertices — which exploits the 
knowledge of densities and the absence of renormalization of the involved amplitudes — 
is the technique of line- dressing, that we describe schematically below. 

Let us suppose that the amplitude A& corresponding to a connected OPI diagram 
D' with V loops is analytically known, with the condition that the mass m'j associated 
to the j-th propagator is different from all other masses in residual propagators, 

r L ' 

A&:= / ]\d% [block] 

let us suppose in addition that the amplitude A^ 2 ' G \p 2 ) associated to some two-point 
connected OPI diagram D admits a spectral representation as the one given in (|45|). 
Then the amplitude of the diagram Q" obtained by replacing (or, familiarly, "dressing") 
the j-th line of diagram D' with diagram D, 

v 

A&, : = / Y[d% [block] Av^ikf), (47) 
^ i'=i 

can be written as 

A v" = f Q dM 2 Pv(M) (Azy\ m >= M ). (48) 

The term [block] stands for the same ^independent expression in both (J46|) and (J47)) . 

Please bear in mind that we used the — implicit but essential — hypothesis that D, 
D' and D" need not to be renormalized. Examples of A$y eligible to line-dressing are 
all the one-loop functions whose expressions are known for non-equal masses (see |41| . 
[2*U]). in all kinematical configurations (see [4"2*]). 

However, it must be noticed that the practical implementation of (J48|) can be very 
complicated; thus, we decided for the moment to make use of the technique of line- 
dressing only in the case of the first two- loop diagram on the left of figure H3 which can 
be seen as a triangle with one line dressed by a bubble: for this effective vertex, the 
resulting cost turns out to be C — 1. 



(46) 



3.4-3. Low-cost subdiagrams In some cases, we can take advantage of particular 
"hidden" symmetries of the diagram we are dealing with to obtain low-cost integral 
representations for some higher-loop functions, with full dependence on external 
momenta. 

An example of such a situation is given by the "square- with- diagonal" diagram of 
figure 121 which — by means of a clever assignement of loop momenta in the spirit of 
section FQ1 — can be parametrized in terms of a one- loop triangular effective vertex and 
only 2 residual integrations. In fact, using the momentum assignement described in 
figure El and performing the integral over £ in spherical coordinates (\£\,cos8,<f)) with 
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Figure 2. "Square- with-diagonal" diagram, and — on the right — the momentum 
assignement and the choice of a triangle effective vertex leading to its low-cost 
parametrization. Please remark that vectors e.% denote the incoming external momenta. 

the z-axis chosen to be parallel to the external momentum e\, it turns out that the 
triangular effective vertex is ^-independent; consequently, the integration over of the 
two residual propagators factors and can be performed analytically. 

The underlying symmetry leading to angular independence is still present in the 
more general situations of diagrams with two or three external legs, respectively, 
bordered by an additional chain of propagators connecting two extremal vertices - 
that is, vertices connected to external lines — ; in this sense, in fact, the "square- 
with-diagonal" diagram can be considered as a triangle bordered by a chain of two 
propagators. Using in such two general parametrization analogous to that in 

figure El and spherical coordinates oriented as explained above, one readily obtains that 
the three-leg (resp. two-leg) subblock is ^-independent (resp. (9, 0)-independent) and 
that these angular integrations involve only the propagators belonging to the bordering 
chain: consequently, such integrations may possibly be done analytically. (As usual, the 
absence of renormalization is an essential hypothesis for this property to safely hold.) 
In the case of a two-leg subblock bordered by a chain the considerations above imply 
that if the angular integral over (8, <fi) of the chain of propagators is analytically known 
one can spare two integrations. The same result can be obtained by the technique 
of line-dressing one-loop skeletons with the spectral density of the two-leg subblock 
in question (see section I3.4.2|) . if this density is available: this latter technique seems 
more adequate for the implementation of two-loop diagrams corresponding to a one-loop 
bubble bordered by a chain of propagators. 

The parametrization described above turns out to be particularly useful when 
applied to the family of two-loop diagrams constructed by bordering with a chain two 
extrema of a triangle one-loop subdiagram; all such diagrams can be evaluated at cost 
C = 2 by performing the integration over <fi of the corresponding chain of propagators, 
and using the known analytic expression for the triangle subdiagram. The simplest 
graphs belonging to this family are shown in figure El together with some useful bordered 
bubbles: employed as effective vertices they will prove important, for example, to reduce 
the complexity of some 7-loop diagrams back to a more manageable size. Following 
considerations exposed above, we will assign a cost C = 1 only to the first graph, which 
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corresponds to a bubble bordered by a two-line chain or, equivalently, to a line-dressing 
of a triangle with a bubble density — see section 13.4.21 — ; even if similar considerations 
could be applied to the second graph, which is as well a bordered bubble, we will content 
ourselves to treat it here as an effective vertex of cost C = 2, since the development of 
the code corresponding to the C = 1 version would require a much larger programming 
effort. 




Figure 3. Principal two-loop low-cost diagrams (no more than two integrations are 
required for their evaluation, see discussion in section 13. 4. 3f) . 

3.4-4- Full-cost subdiagrams Quite surprisingly, implementing effective vertices as 
standalone basic blocks can lead to important simplifications even if we are not able 
to compute the effective vertex in a cheap way and spare integrations; the mere fact 
that more than one basic block can be replaced in one amplitude is sometimes enough 
to give rise to reductions in its final cost. 




(a) (b) (c) (d) 

Figure 4. A difficult graph together with its treatment. 



An enlightening example of such a situation is offered by the diagram shown in 
figure |3J a; due to its complexity it looks formidable, until one realizes that it can be 
split as in figure 0]b and expressed in terms of two identical blocks, each having the 
form illustrated in figure EJc; thus, the block of figure EJc can be considered as a new 
kind of effective vertex with 5 asymmetric legs, which can be identified and replaced 
into larger diagrams in the usual way (as illustrated in figure Eld for the case of diagram 
of figure |U a itself). 

From a more extensive analysis of figure EJd — using the definition of cost (J54j) to 
deal with the nonzero-cost effective vertices, and computing loop costs as described in 
section 0] — we obtain for this diagram a final cost of 6 (3 for the computation of loop 
integrals, plus 3 for the computation of the effective vertices, which is in fact the full 
cost for computing such two- loop diagrams); on the other hand, it should be noticed 
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that due to topological obstructions only two one-loop ordinary "low-cost" functions 
could be replaced into the diagram of figure E|a if we did not introduce the "full-cost" 
effective vertex of figure EJc, leading to a much larger final cost of 9. 

3.5. Momentum representation: identities for vanishing external momenta 

In this section we collect some simplifying diagrammatic identities, which take advantage 
of the fact that our amplitudes are always computed at zero external momenta (and 
that renormalization is trivial in our scheme). 

3.5.1. Factorizable amplitudes Quite often connected OPI diagrams in (0 2 ) 2 theory 
show cut-vertices, i.e. — in this context — vertices such that the separation of their four 
legs into two appropriate groups of two legs disconnects the diagram in two parts. In 

Figure 5. Factorization property of amplitudes of diagrams with cut- vertices when 
external momenta vanish. The vector refers to the vanishing of the momentum flow 
between the two blocks C\ , C%. 

figure El we depicted a cut- vertex v, which divides its graph in two subgraphs C\ and C 2 . 
At vanishing external momenta, momentum conservation can be obeyed if and only if 
there is no total momentum flow from C% to Ci] using the locality property of momentum 
representation and an appropriate choice of loop momenta we can always separate the 
integrations relative to block C\ from those relative to block C2, expressing the original 
integral as the product of two integrals, each one associated to a lower order connected 
OPI diagram evaluated at zero external momenta and carrying a J <p 2 insertion in place 
of the original (</> 2 ) 2 vertex v. In the following we will refer to diagrams with one or 
more cut- vertices as to factorizable diagrams. 

It turns out that a relevant fraction of the diagrams needed to compute RG- 
functions — about the 20% up to the 7-loop level — is factorizable in terms of lower-order 
diagrams; these diagrams are in turn either easier to evaluate than the original one, or 
already known. This idea is fully implemented in our framework (the corresponding 
logical block has been labeled as "filter for chain diagrams" in figure EH]). 

3.5.2. "Quail's leap" Other interesting relations can be obtained at zero external 
momenta when additional symmetries are present in the expression of the amplitude 
for a given diagram. In fact, it is sometimes possible to look at the diagram in terms of 
a composition of blocks of propagators, with a momentum flow such that blocks can be 
rearranged; one can in this way obtain other diagrams, which differ in structure from 
the starting one but whose associated integrals have the same numerical value at zero 
external momenta. 
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Figure 6. An example of "quail's leap" identity. The dashed lines express here 
momentum flow between blocks, and do not correspond to real propagators. 



In figureElwe give a pictorial illustration of an example of such set of identities; from 
now on, we will familiarly call it the "quail's leap". The analytic counterpart of figure |H1 
follows straightforwardly if one writes down the explicit form of amplitudes A® and A& 
as schematically defined by figure EH and then relates them using the commutativity of 
the ordinary product in momentum space; as stated, this leads to the identity 



(It should be noticed here the essential role played by the vanishing of external 
momenta; this condition automatically enforces momentum conservation in the diagram 
D' associated to the amplitude with permuted blocks.) In spite of the simplicity of its 
analytic derivation, consequences of "quail's leap" are non-trivial at all. Starting from 
the identifications 



we can easily show an example of use of the relation in figurelBl substituting the indicated 
blocks into it we readily get to the following non-trivial diagrammatic identity 



Strictly speaking, it is not guaranteed that the diagrams obtained applying identities 
such as the one illustrated in figure will belong to the same correlation function as 
the diagram one starts with; in practice, it turns out that the number of such useful 
cases does not represent a large fraction of the diagrams required, and this trick is not 
implemented in the current form of our framework. 

3.6. Other aspects of momentum representation 







In this section we collect a couple of additional techniques, which, although not 
leading to a direct reduction in the dimensionality of the integrals to be performed, 



Automatic parametrization of amplitudes 



25 



are nonetheless essential to complete our task of parametrizing amplitudes. In both 
cases, anyway, we take advantage of momentum representation of amplitudes to put 
the required manipulations in a neat diagrammatic form. As usual, triviality of 
renormalization in the I-scheme is tacitly exploited. 



3.6.1. The derivative with respect to p 2 In this section we show how to compute 



<9if' 0) /dp s 



p=0 



from the diagrammatic expansion of If 



(2,0) 



for a better insight, we 



keep a generic dimension d in the algebraic expressions presented below. 

The derivative w.r.t. p 2 of a function f(p 2 ) can be expressed in terms of derivatives 
w.r.t. momentum components p M by using the chain rule. In particular, the following 
relation 

d d „ , ox _ . d „ , on . o d 2 , , 2 , 

:/ (P ) 



— f ( p 2\ = id— f (p 2 ) + Ap 2 - 

dpa dp a v ' dp 2 v ' dp 2 dp 2 ' 



turns out to be useful. Assuming in addition f(p 2 ) to be smooth at p 2 = 0, it follows: 



1 

2d 



d d 



p=0 — \dPnd Pfl 
The derivative w.r.t. p 2 of rp' "* 



(50) 



p=0 



at p = is obtained by distributing the linear 
differential operator d/dp 2 over each Feynman amplitude contributing to r| 2 '°' ) at 
nonzero p, and then applying the identity (|50|) to each amplitude. (It should be noticed 
that each amplitude is a function of p 2 only, and is regular at p = because all masses 
in propagators are nonvanishing.) Commuting the derivatives with loop integrals leads 
us to the problem of evaluating the action of differential operator d/dp^ -d/dp^l 
on a product of p-dependent propagators (obviously, p-independent propagators do not 
take part in the process, and stay unchanged). In the following we will stick to the 
practically convenient case when the external momentum p flows without fractioning 
from one external point of the diagram to the other through a chain of propagators (the 
extension of given formulas to the general case is straightforward). 

We label with the index a the lines belonging to the chain, and with the index b 
the lines not belonging to it; calling k a and fc& the internal momenta associated to each 
line, and using for shortness the notations Pr(fc) := 1/(1 + k 2 ), Pi^(k) := k^/(l + k 2 ) 2 , 
we easily obtain: 



l__d d_ 

2d dp. dp h 



Y[Pi(p + k a ) l[Pr(k b 



Y[Pr{k 



+ 




p=0 

^ Pr(k a ) 2 - - d Pv{k a f 



Pi^k a )P^{k 



n pr ^« 



(51) 



n pr (^ 

Remarkably, due to the freedom we have in choosing how to parametrize amplitudes the 
final integrated result will not depend on the choice of the chain; however, it is evident 
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from (J51|) that to compute the derivative of an amplitude using a chain of length / 
one must deal with I (I — l)/2 + 2Z generalized diagrams — involving, to worsen things, 
additional two-leg vertices or vector propagators. As a consequence, to get an efficient 
automated framework for parametrizing amplitudes one must add a stage where all 
diagrams of rp' ^ having a nontrivial dependence on external momentum are scanned; 
since each chain allows the computation of the derivative w.r.t. p 2 of the original diagram 
via the values of the set of graphs induced by ()51|) . the goal will be to find the chain 
which is "minimal" in terms of the set of integration costs of such induced set of graphs. 

This processing stage has been indeed implemented in our code, and is in fact 
responsible for the majority of the time spent to parametrize the 2-point function: the 
large number of chains to be examined, together with the large number of diagrams to 
be parametrized for each chain, give in general rise to an unpleasant explosion of the 
number of diagrams to be examined. 

We notice as a final remark that a special prescription can be formulated to compute 
the derivative of factorizable amplitudes which we examined in section IB. 5. 11 in fact, in 
this case we can arbitrarily decide to confine the chain of propagators intervening in (Jo"Tj) 
to the first block of the diagram, that is, to the factored block to which the two external 
legs of the diagram are connected; in this case, we will obtain the derivative of the 
whole diagram by computing the derivative (hopefully simpler) of the first block, and 
then multiplying each of the terms obtained from the r.h.s. of (151)1 by all the remaining 
original factored blocks — which will be left untouched by the action of the differential 
operator d/dp 2 , since due to our particular choice of the chain they do not contain any 
of the propagators involved in the derivation. This consideration explains the intricated 
structure of the upper part of figure 13.11 

3.6.2. Insertions of J <p 2 The diagrammatic expansion for r : (2,1) (p = 0,q = 0) can be 
obtained from that of rj 2 ' '* (p = 0) by inserting J into each graph which belongs to 
the 2-point function; from a diagrammatic point of view this operation corresponds to 
the insertion of a two-legs vertex (with no additional external momentum), and can be 
performed by replacing in all possible ways one propagator of the original diagram (that 
is, one line) with a chain of two propagators; when the same final graph is generated in 
many different ways, its multiplicity has to be taken into account (see figure [7| for an 
example). 




Figure 7. A an example of diagram contributing to T\ ' , and the corresponding 
groups of diagrams in r[ 2 ' 1 ' > generated by the insertion of one J ^ operator. 



As observed in [221, it turns out from simple topological considerations that these 
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diagrams correspond to a subset of the diagrams of If (p = 0); this fact implies 
that we do not need to recompute the numerical values for the integrals associated to 
diagrams of Tp '[P = 0,q = 0) at some loop order L, since such values can be deduced 
from a lookup in the list of values already computed during the evaluation of the 4-point 
function at L loops — this nice property cannot be extended, however, to corresponding 
symmetry and O(N) factors. 

4. Automatic parametrization of amplitudes 

We are now ready to examine in full detail how an automatic framework for finding out 
minimal-cost parametrizations can be built. 

4-1- Overview 

Provided that the key concept of the method is the substitution of known effective 
vertices into the analytic expression of the amplitudes, the problem is how to fully 
exploit such an idea in a systematic, automated and optimal way. 

The starting point will be a choice of a set E := {Si, . . . , £jv cv } of effective vertex 
functions. To each effective vertex S e is associated a positive integer C e called cost. When 
S e is computed by use of an appropriate integral representation the cost is defined to 
be the dimensionality of the integration, see ()53|) ; otherwise, if the vertex function is 
known analytically in terms of elementary functions, the cost C e is zero. 

Formalizing now the ideas of section 13.41 once N ev effective vertices have been 
replaced into the amplitude A® corresponding to a given Feynman diagram D, with 
L res residual loops remaining, the new expression for that amplitude will be given by 

/Lrcs / Ires -i \ / ^ev \ 

(Jii+^J (n 5e(sPe) J (52) 

where the effective vertices with non-zero cost are expressed as: 

£ e (§P e ) := J d% S e (<SF e ,Q . (53) 

7 res above is the number of residual propagators and kj the momentum flowing in each 
of them; the set SP e C { l^.l^ \ li,l 2 G 1 . . . L Tes with h < l 2 } encodes the dependence 
of the effective vertex S e in terms of scalar products among loop momenta and plays 
a key role in the choice of the most appropriate angular measure, as will be shown in 
section 14.41 

The inclusion of effective vertices with non-zero cost in the set of possible 
replacements imposes a refinement of the definition of the cost function itself. The 
problem in question is how to extend the definition of total cost as dimensionality of the 
involved integration — which is appropriate when the amplitude is expressed in terms 
of one single multidimensional integral — to situations in which two or more effective 
vertices with non-zero cost — and, thus, two or more blocks known only in terms of some 
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integral representation — are present in amplitude parametrization and are themselves 
integrated as in equation (|52[): the solution to this problem is somewhat implementation- 
dependent. In the following we will refer to the concrete situation in which the numerical 
integration of the parametric integrals defining the blocks with non-zero cost via (|53|) 
is performed separately for each block, and also separately from that of loop momenta. 
In the presence of two or more blocks with non-zero cost it seems quite natural - 
and convenient — to prefer this type of implementation to that based on the opposite 
strategy of collecting and numerically evaluating all integrations together in a single 
multidimensional integral. 

Following these lines, we will assign to an amplitude parametrization like the one 
shown in (J52|l - (J53*Jl the refined total cost 

C := ( loop cost ) + ( max C e 

V / \l<e<N cv 

The "loop cost" appearing in this expression is given by the still unknown residual 
dimensionality of the integration of loop variables Yii=i d 3 £i- It will be the goal of 
section POl to compute, and possibly minimize, this number. 

Some additional comments about formula (J51j) are in order. First of all, this refined 
definition of cost can be easily justified by trusting equation (J32|) and supposing that the 
number of operations of a sequence of integrations is dominated by the integration(s) 
of highest dimensionality. Secondly, a finer structure exists: in essence, when a 
parametrization has a total cost C > we can still say that to evaluate it we roughly 
need a number of operations proportional to those required by a C-dimensional integral; 
however, 

• if we adopt definition (JH%|) for the total cost, a constant proportionality factor in 
the number of operations is missed in all cases when a parametrization requires the 
evaluation of more than one single effective vertex having a cost exactly equal to 

maxi< e <jv ev C e 

• even in the simple case when only one effective vertex needs to be integrated, 
and we are then free to choose whether or not to carry out the integration of the 
effective vertex along with that of loop variables, from the point of view of numerical 
integration the task of computing /(...)• J„ £ e is not entirely equivalent to that 
of computing f Cx ~ (...) £ e (S being the integration domain of the effective vertex 
with non-zero cost, and £ the integration domain of loop variables); this happens 
first of all because more efficient non-separable integration rules can exist allowing 
to integrate the domain £ whole, and then because a higher intermediate 
precision is in general needed for the evaluations of J" (. . .) and f~£ e if we want to 
obtain a given numerical precision for / £x= (• • •) £ e - 

Nonetheless, since the task of specifying a precise formula for the computation of the cost 
is indeed a very difficult one — in the light of considerations carried out in section 13. 31 
we will stick to simpler formula (JHlj) throughout all this article. 
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Consequently, having picked up a set of known effective vertices E with associated 
costs {Ci, . . . , Cjv ev }, our strategy to reduce the total cost C for a given diagram D will 
be as follows. 

We start by considering the associated amplitude A® parametrized in momentum 
space, 

j 1=1 »=i 1 + Ki 

Then, we try to reduce the cost (p>4"|) by applying the following guidelines: 

(i) we look among all possible maximal substitutions in the graph D of subdiagrams 
from the given set of effective vertices E. Each replacement is encoded via an 
effective graph D*, which contains as many new vertices as possible describing the 
replaced effective vertices, and gives a new expression for the amplitude in terms 
of fewer integrals. An algorithm to scan the space of all possible substitutions is 
presented in section l4~2l 

(ii) we reduce the number of scalar products ti x .t\^ our integrand function depends on 
by means of a careful choice of loop basis (see section 14 .3|) 

(iii) we maximize the number of trivial angular integrations by selecting the cheapest 
angular integration measure which is allowed by our previous choices of substitutions 
and of loop basis (as explained in section 14. 4 j) 

(iv) while picking up the best angular measure, we maximize at the same time the 
possibilities of factoring away simpler integrals, like for example 

/Xttf) (lT(^)"> 

(see as before section [Ol for more details). 

It should be noticed that different (and perhaps more general) strategies can in principle 
be devised; of course, the list of tricks to minimize the number of residual integrations 
we just presented is only a faint approximation of what a well-trained "human neural 
network" can do by examining all diagrams one after the other, and finding out what 
is likely to be the real optimal parametrization. However, since teaching mathematical 
intuition to a computer is clearly impossible, what we are interested in here is the 
much simpler task of creating a fast and robust algorithm relying on a reasonably small 
number of simplifications, and working reasonably well for the majority of diagrams. A 
more in-depth discussion about this philosophy and the optimality of the solutions it 
gives can be found in section El 

4-2. Substitution of effective vertices 

To fulfill our purpose of obtaining parametrizations of a given diagram which have a 
minimal cost in the sense of equation ([54)1 . we must as a first step find out all possible 
maximal substitutions of effective vertices for that diagram. 
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Algorithm 1 (Tree of Substitutions) To find out all possible maximal substitutions 
for a given Feynman diagram D, 

(i) read the set E of known effective vertices, graphically expressed as pairs 
p e = { subgraph s e , effective graph e e corresponding to replacement } 

(ii) define 

(a) the tree of substitutions: to each node is associated a diagram which is obtained 
by replacing a collection of subdiagrams of D — non necessarily maximal— 
with corresponding elements of {e e }; the node also contains the list M. of 
all remaining possible replacements, each one accompanied by a status flag 
describing whether the related alternative has already been visited or not 

(b) the set M of maximal substitutions found so far 

(c) the lookup node 9t* 

(Hi) initialize: 

(a) create the root node yi of the tree of substitutions; associate to it diagram D, 
and a list M(0to) appropriately computed 

(b) m < — 0, m* < — m 

(iv) if the next unvisited replacement R in the list of replacements R(W) exists, 

(a) updateR(W) 

(b) apply R to the diagram associated to lookup node 71*, thus obtaining a new 
diagram D' 

(c) attach a new node 71' to 71*; associate to it diagram D' , and a list IR(W) 
appropriately computed 

(d) update the lookup node: 71* < — TV 

(e) go to step 
else 

(a) put the diagram associated to lookup node 71* in the set M of maximal 
substitutions 

(b) do backtrack up to the nearest node whose list of unvisited replacements is 
not empty; 

if such a node exists, 

1. set that node as the current lookup node 

2. go to step Hv]) 
else end. 

As an example, we now illustrate how algorithm 1 operates when it is applied to 
the diagram D and the set of effective vertices E which are presented in figure |H1 

In figure IH1 we can observe how the traversal works: the algorithm always tries to 
push up the level of the tree, in the hope of finding out one effective vertex more in 
the same lookup diagram. A view of the entire tree (more sketchy) can be found in 
figure [TUl It turns out that for these choices of 2) and E there are exactly four possible 
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S := 




E 




Figure 8. A test case for algorithm Q: a diagram D and a set of effective vertices 
E to be replaced in it. 




. [0][50][. 



Bubbles 

[1][49] 



Triangles 
j511[...1 



Squares Pentagons 



[2][4S] 




[52][...] 



[10] [12] [14] [16] [18: Accept] [22][24] [26] [28] [30: Accept] [34] [36] [38] [40][42: Accept] 



3u Tr Sq Pe 
[11] [13] [15] [17] 



V V 

Bu Tr Sq Pe 

[23] [25] [27] [29] 



7 V 

Bu Tr Sq Pe 
[35] [37] [39] [41] 



Figure 9. A partial tree of substitutions for the example of figure |HI The order 
of traversal of each node is explicitly specified between square parentheses. A cross 
denotes the fact that the corresponding substitution is not possible. The label "Accept" 
means that the set of replacements associated to the node is maximal. 



different maximal substitutions, three with a residual loop number L res = 2 and one 
with L res = 3. 

Some general remarks are in order here: 

(i) we can easily convince ourselves that due to possible topological obstructions 
the distance from the root of each terminal leaf of the tree of substitutions is 
not the same. From a different perspective, we can say that the final maximal 
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Figure 10. The complete tree of substitutions for the example of figure [S] The 
numbers in frame corners indicate the element of the basis of effective vertices which 
is being replaced into the diagram in each case. 



substitutions will have in general a different number of residual loops, thus leading 
to parametrizations with a different number of final integrations. The importance 
of such a conclusion in the perspective of the calculation of cost C should not be 
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underestimated 

(ii) by taking a closer look to figure [TUJ it turns out that the tree built up by this 
naive formulation of the algorithm must be pruned in order to be used in effective 
calculations. In fact, it is clear that the same final parametrization can be in general 
reached following a lot of different branches of the tree: being subdiagrams with 
a given form but with a different embedding in the diagram distinguishable from 
each other, the same final set of n vertex functions can be obtained from exactly 
n\ different tree paths. Adding to this problem the observation that this algorithm 
too is plagued by the problem of isomorphic copies, one can very well understand 
how the size of such a naive tree is soon pushed to astronomical values. Luckily, the 
use of ordered searches for vertex functions brings the complexity of the algorithm 
back to a manageable size. 

4-3. Choice of loop momenta 

After a choice of a maximal substitution of effective vertices in the original diagram 
D has been performed, resulting in an effective diagram D*, the second step in the 
generation of parametrizations is clearly the task of assigning loop momenta to D*. 
In spite of the fact that the choice of loop momenta is often looked at as a trivial 
operation, a favourable assignement of loop momenta — if complemented by successive 
optimization of angular measure and maximization of simple angular integrals over 
residual propagators, see section 14.41 and section 14.4.31 — can in many cases provide 
further reductions in the number of final residual integrations. 

The following facts — which are common lore in graph theory — will be very useful to 
get to algorithmic which allows the sequential generation of all possible loop assignments 
(further references and definitions adapted to more general contexts can be found for 
example in [H], j3S] and |36j). 

Definition 1 Given two subdiagrams S\ and S2 of a diagram D, we define s = Si © S2 
as the subdiagram of® whose lines are present in s% or in s 2; but not in both (exclusive 
sum of lines). 

Definition 2 A ri-cycle — or n-loop in physicists' terminology — is the diagram 
composed by the set of vertices {t>i, . . . , v n } and the set of distinct lines 
{v\ — V2, ... , v n — Vi}. A tree is a connected diagram with no cycles. A spanning 
tree for a connected diagram is a tree subdiagram with the same vertices as the original 
diagram; by topological count, the operation of adding a single line to a spanning tree 
does always produce a subdiagram with loop number equal to 1. 

Theorem 1 Given a connected diagram D, its loops span a vector space Yl in Z2 w.r.t. 
the operation ©. The dimension of such a vector space is given by L, the number of 
loops of the diagram familiar to physicists. 

Not all the elements of Yl are loops, but all the loops of the diagram are contained 
in Yl] thus, this theorem gives us a powerful tool to enumerate them all. 
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Theorem 2 To get a loop basis B for the vector space Yl associated to a connected 
diagram D, 

(i) produce a spanning tree X for the diagram 

(ii) for each remaining line I G D such that t G 3) but 1 G" X, consider the subdiagram X[ 
obtained by adding the line t to 1, that is %i = TUl; produce a cycle q by identifying 
it as a subset of%\. 

The set B = {ci, . . . , Cl} of cycles generated in this way is the desired basis. 

As a consequence, we can now formulate the following 

Algorithm 2 (Loop bases). To produce all possible assignments of loop momenta in 
a given connected diagram D 

(i) build a spanning tree X for D 

(ii) deduce a first loop basis Mi for the loops in the diagram 2) by applying theorem^ 
to the couple {D, X}/ the property L = card(Bi) will hold 

(Hi) deduce the set L of all loops contained in the diagram by calculating all possible 
linear combinations over Z 2 of the elements of B x and by discarding combinations 
that are not loops 

(iv) for each extraction of L loops c^, . . . , Ci L G L 
if loops Cjj , . . . , Ci L are linearly independent 

(a) take the set {c ii; . . . , c iL } as a new possible choice of loop momenta 

(b) for each c in 

1. assign an arbitrary orientation to it 

2. increment (according to the chosen loop orientation) the momenta 
associated to the lines of Ci n by a quantity a n £ n - - where the a n ; s 
are arbitrary non-null coefficients. 

As an example of application of algorithm El we illustrate in detail the loop 
parametrization of the "cat's eye" diagram; it is drawn in figure ITT1 together with a 
choice of a spanning tree, T, and with the loop basis Bi obtained by applying theorem[2] 
to T. 

The three cycles which are elements of B x are identified (from left to right) with the 
following basis vectors: {0,0,1}, {0,1,0} and {1,0,0}. The set of all cycles contained 
in the original diagram is obtained by working out all possible linear combinations of 
such basis vectors in and by discarding the resulting graphs which have more than, 
or less than, 1 loop (see figure IT^|) . All possible loop bases will then be given by all 
choices of three linearly independent loops among the six so far produced: for instance, 
B 2 = {ci, C2, C5} is a valid basis but the same is not true for the choice {ci, c 2 , C3}, due 
to the fact that c 3 = Ci © c 2 . 

Then, given a choice of loop basis and of suitable orientations, momentum 
assignements for the lines will follow by attributing momenta otii-i, a 2 £ 2 and 03^3 (with 
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Figure 11. "Cat's eye" diagram: a spanning tree T and the associated loop basis Bi. 



{0,0,0} Cl = {0,0,1} c 2 = {0,1,0} c 3 = {0,1,1} 



< 1 loops 




c 4 = {1,0,0} C 5 = {1,0,1} Ce = {1,1,0} {1,1,1} 




> 1 loops 



Figure 12. "Cat's eye" diagram: all cycles. 




Figure 13. "Cat's eye" diagram: loop assignment deriving from loop basis B 2 = 
{ci,c 2 ,c 5 }. 

an 7^ arbitrary) to the three cycles of the basis, respectively: in figure [T3] we depict 
the loop choice which is obtained for the "cat's eye" by choosing the B2 defined above. 

It should be noticed that even the computation of loop bases can give rise to a 
combinatorial explosion of generated data, in particular when the number of residual 
loops L res of the effective diagram D* is large: this fact can be easily understood by 
realizing that each loop basis is an extraction of L res linearly independent elements out 
of a set L composed by up to 2 Lrcs — 1 elements. 

Finally, a summary of the entire parametrization procedure for the diagram 
introduced in figure |H1 is presented in figure H3| Please remark that the deduction 
of the costs shown in this figure relies on the choice of angular measures, which is the 
subject of the next section IPl 

4-4- Choice of angular measure 

In this section we deal with the problem of explicitly parametrizing the integration 
measure for an integral expression of the form 
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Choice of angular measure 



Figure 14. Example of figure[S]after various choices of loop momenta. The costs will 
be deduced only at a later stage, after angular measures for each diagram have been 
picked up (see section 



where, using notation previously introduced in section 14. 1[ F depends" 1 " on a subset 

SP of all scalar products { t\ x .£i % \ h,h £ 1 • • • -^res with h < h } among integration 

+ In reality, the fact that the L rcs x L rcs Gram matrix of scalar products (ii .£j ) has rank not greater than 
the integer space-time dimension d implies — when L rcs > d — the existence of (L rcs — d)(L res — d+l)/2 
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vectors. 

Since an hyperrectangular domain seems to be the most convenient set-up when 
we are to perform multidimensional numerical integration, in the following we will 
restrict ourselves to the particular choice of spherical coordinates, and will accordingly 
decompose loop vectors as the product of their norm times a versor: t\ = l\ t\. 
Consequently, given an integral 

Y[d£i l\ dtti F^SP) , (55) 
1=1 

where 

SIP C { \2 h =: cose hl2 \l h l 2 el... L vcs with h < l 2 } (56) 
l=(*i, •••^0> (57) 
we examine here the problem of explicitly writing down an angular measure Yli d£li 
which is minimal in the number of required residual integrations, that is, which allows 
a maximal number of angular integrations to be carried out analytically. Although the 
spirit of this idea can be easily generalized to a different number of dimensions d, our 
explicit formulas will be specialized to the case d = 3. 

This investigation is motivated by equation (|52|) . which shows that our most 
general amplitude is precisely of the form just mentioned; the cost (jB^) associated to a 
parametrization of a given diagram cannot be fully evaluated without knowing the cost 
of the integration over loop momenta, and due to this fact one is forced to explicitly 
determine such a "minimal" angular measure. 

It should be noticed that parametrizations of angular measures completely different 
from the spherical ones can of course be considered — we think e.g. to the techniques 
used in [46J — , possibly leading as well to a significant reduction in the number of 
residual integrations; however, we have not considered them in this work, due mainly to 
the fact that high-performance numerical integrators seem to get along best with simple 
hyperrectangular domains. 

4-4- 1- Cosine diagrams As we will see later on, §P, the subset of scalar products 
on which function F in equation ()55|) depends, plays an essential role in determining 
the optimal choice of angular integration measure. This information is encoded in the 
so-called cosine diagram associated to function F, (£, defined as follows: 

functional relations among such scalar products, and consequently that at most dL Tes — d(d — l)/2 
scalar products out of L Tes (L Tes + l)/2 can be functionally independent. As a consequence, the same 
function F might be defined in terms of different expressions, each involving an a priori different subset 
of scalar products (li.lj). For shortness, in the following we use the words "a function F" meaning 
"an expression — fixed throughout the analysis — which gives a realization of the function _F" ; however, 
notice that the minimality of the set of scalar products on which the chosen representation depends, 
SP, is never assumed in our approach. In particular, eventual functional dependencies among the scalar 
products are automatically taken into account in our framework, being the parametrizations of scalar 
products obtained from parametrizations of vectors — see eq. (|58J) . 
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(i) for each integration variable £[, add a vertex v\ 

(ii) if function F carries an explicit dependence on ii x .li 2 (with li < / 2 ), draw a line 
connecting vertex vi ± to vertex v\ 2 . 

At this level the knowledge of SP is equivalent to that of (£, but further information - 
encoding additional properties of F — can be attached to the cosine diagram by allowing 
different kinds of lines in (£. A fruitful application of this technique to the specific class of 
function we are interested in is postponed to section 14.4.31 In the following section 14.4.21 
we will discuss at a more general level the problem of a "good choice" of parametrization 
of the angular measure. 



4-4-2. Choosing axes The abitual way of parametrizing the angular measure Yl l dfli is 
to decompose each dfli as 

dfli = — d(cos6i) dtfii] 

however, this decomposition requires a previous choice of an axis z\ — in terms of which 
angle Q\ is measured — and of an axis x\ — such that angle (pi can be defined as in 
usual 3-dimensional spherical coordinates. (Please remark that in the following xi, z.\ 
will denote a choice of non-degenerate — but not necessarily orthogonal — axis versors.) 
Once a choice of axes has been picked up, one can explicitly parametrize all the domain 
of variation of versor £i in the standard way as 

li = cos 9\Z\ + sin 9 t (cos (pi (z t A jfcj) Azj + sin 0/ z x A Xij (58) 

with 9i e [0, it) and (pi E [0, 2tt) 

(our standard notation v := v/|v| has been used here). In addition, after a 
parametrization is given for all versors ii one can usefully encode this information in 
the matrix of scalar products (MSP), which is defined as 

, v f 4 2 h (e, (p) if i< h < i 2 < L res , 

(MSP) W2 (9 A) := Zl /2U '^ - 1 2 - res (59) 

10 if 1 < <2 < tl < -t^res- 

(Please remark that 9 means "an appropriate subset of 9 1 , . . . , 9 Lvcs " , and the same 
holds for (p; refer to Appendix A for an explanation of our notations.) 

Now, supposing we have formulated a choice of axes Xi, zi for each t\ and prepared 
the matrix of scalar products (MSP) associated to this parametrization, the question 
arises: how many nontrivial angular integrations are necessary to parametrize the 
angular measure Yli d£li for a given function F which depends on a (sub-) set of scalar 
products §P via a specific cosine diagram £? 

In present framework the answer follows straightforwardly: just use the 
informations encoded in the cosine diagram € and in the matrix (MSP) of parametrized 
scalar products to build up the set AV containing the "nontrivial" residual angular 
variables which F really depends on: 

AV := [J {9,(p\ {MSP) hh depends on variables 9,$}. (60) 

line h—fa £<£ 
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The relevance of the set AV is immediately apparent: given an angular measure J7, dVti 
induced by a specific choice for a parametrization, all angular variables 0, involved in 
its realization which are not elements of AV correspond to flat directions, and can in 
accordance be trivially integrated out. 

Coming back to the issue of parametrization via (|58|) . it clearly appears that a great 
arbitrariness is still left over: in fact, given a function F with cosine diagram £, many 
choices can be imagined for versors Z\ and Xi, each one leading to an a priori different 
number of final angular integrations, depending on the particular form of (£. 

One result of the present work is the classification of all possible cosine diagrams 
relevant for integrations up to L res = 4 and the systematic development of appropriate 
choices of axes to lower as much as possible the number of final angular integrations. In 
the following we will clarify the theoretical framework, and introduce the reader to our 
parametrization strategy by giving some examples. 

Let us suppose of having to cope with two integrand functions, Fi, F 2 depending 
respectively on four loop versors £i, ... , £ 4 via the set of scalar products SPi and SP2 
defined as: 

SPi := I (-1-&2-, ^2-^3, } (61) 

SP 2 := { \%, \% } (62) 



As a warm-up we consider the natural choice of axes parametrization consisting in 
picking up two fixed nondegenerate external vectors Z and X, and in parametrizing all 
the versors i\ w.r.t. these axes via (|58j) : the following table encodes this possibility. 



Axis choices 


% .. 






Z 


z 




X .. 


. X 



(63) 



The choice in (|U3*|) would lead in the L res = 4 case to the particular form of (MSP) given 
by 

( c0ic0 2 + s6»is6» 2 c(0i - 02 ) c9xc9 3 + s6'is6'3c(0i - 3 ) c9ic9 A + s6>is6» 4 c(0i - 4 ) \ 
c9 2 c9 3 + s9 2 s9 3 c((f)2 - 3 ) c9 2 c9 A + s9 2 s9 A c((f) 2 - 4 ) 

C0 3 C#4 + S0 3 s6 , 4 c(0 3 — 4 ) 

VO / 

Following equation ([BOjl we readily obtain the sets AVi 5 fi xc d and AV 2i fi xe d of nontrivial 
angular variables induced by the parametrization (|U3*j) for the two cases F 1 , F 2 : 

AVi jfixcd = { 0i, 2 , 3 , 4 , 4>l, 02, 03, 04 } (65) 
AV 2 , fixcd = { 01, 2 , 03, 01, 02, 03 } • (66) 

According to the above results we are now able to estimate the cost of residual angular 
integrations as induced by the parametrization KUty) to 8 for F\ and to 6 in the case of 
F 2 . Impressively enough, the developed formalism reveals quite immediately that due 
to the simpler angular dependences of F 2 2 angular integrations can be spared w.r.t. the 
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"full cost" case of F\. (In fact, a deeper analysis of (|64j) would show that an additional 
integration could be spared for both F\ and F 2 since (MSP) depends only on differences 
of variables; however, this point is inessential in the present discussion because — as 
shown below — one can do much better than that.) 

The key that opens the door to much more efficient parametrizations is the 
possibility to choose as parametrizing axes z t and X\ for £ t appropriate "moving" 
reference vectors. 

A fruitful example of such a choice is the parametrization defined below, which we 
call "standard": 



Axis choices t\ 


% 


% 








% 


il 


% 


il 


Xl 




% 


% 


% 



(67) 



(68) 



in the case of L res = 4 our "standard" choice leads to a cosine matrix (MSP) of the form 

/0 c9 2 c9 3 c9 4 \ 

c9 2 c9 3 + s9 2 s9 3 c(f)3 c9 2 c9^ + s0 2 s#4C04 

C#3C#4 + S6 , 3S6 , 4 C(03 — 4>4) 

Vo o / 

Please remark that in the standard parametrization the expression of (MSP) in dHHI) is 
left unchanged by any legal choice of Xi, Z\ and x 2 . that is why no explicit choice is 
given for such versors in definition ([67)1 . 

The sets AVi jSta ndard, AV2 )S tandard induced by the standard parametrization (JBTj) for 
our examples F±, F 2 follow readily: 

AVi ,standard { #2, #3, #4, 03, 04 } (69) 
AV 2)Sta ndard = { 9 2 , 9 3 } , (70) 

resulting in 5 residual angular integrations for the full-cost case F\ and — quite 
spectacularly — in only 2 residual angular integrations for the simpler case of F 2 . 

From a closer inspection of ([67)1 . and of its implications on the matrix of scalar 
products (MSP) — think to (jQ5|) and its generalizations — it should appear evident that 
standard parametrization implements in an explicit way the possibility of making use 
of the overall SO (3)-symmetry of the integrand to drop the trivial integration over a 
global SO (3) rotation of the t\. In the full-cost case of a function depending on L Tes 
vectors t\ with complete cosine diagram £, the standard parametrization would give 

AV s t a ndard = { 9 2 , . . . , 9 Llcs , 3 , . . . , Lrcs } , (71) 

i.e. 2L ros — 3 residual angular integrations to which one must add L res integrations over 
the vector norms, obtaining the total cost 3L res — 3f, in agreement with the <i-dimensional 
formula introduced in section 13.31 (By the way, a d- dimensional generalization of (|57jl 

f It is remarkable that 3L res — 3 is equal to the maximal number of functionally independent scalar 
products which can be obtained from L rcs 3-dimensional vectors, see footnote at the beginning of 
section POl 
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can be straightforwardly derived. See jUj for a similar result from a quite different 
parametrization. ) 

The example of function F2 above teaches us that when the angular dependence 
of the function is simple, additional integrations can be spared — within the framework 
of standard parametrization — compared to the worst-case equation (J7T|) : however, to 
reach this goal, the properly permuted version of the standard parametization (joTj) must 
be used, think e.g. to a function F 3 with SP 3 := j £1.^3, £2-^3 j- 

All in all standard parametrization (|fi7|) proves to be useful to treat a large variety of 
angular dependencies; however, it should be emphasized that considering other choices of 
moving axes xi, z\ for specific cosine diagrams can sometimes result in parametrizations 
even cheaper than the standard one. 

Let us examine, for instance, the example of an integrand function F4 depending 
on £1, ... , £ 4 via a set of scalar products : = j t\.li, £i-£s, ti.t± j. In this case 
the standard parametrization (|67j) would give a set of residual variables 

AV4 )Sta ndard = { #2, #3, #4, 04 } , (72) 

corresponding to 4 angular integrations; this result could not be improved by 
permutations of (jHTjl . 

On the other hand, the "non-standard" choice 



Axis choices t\ 


% 


% 






% 




% 


Xl 




% 


% 



leads us to the minimal set of angular variables 

AV4 inon -standard = { #12, #13 5 #24 } , (73) 

corresponding to an angular cost of 3: two angular integrations are spared w.r.t. the 
standard parametrization (please remark that in (jTHj) we defined as usual 9^i 2 : = 
arccos ^i^.i^j). 

The variety of parametrizations based on moving axes is enormous, and limited only 
by the imagination of the user. The only important constraint when giving a choice of 
axes Xi, Zi in terms of integration variables £i is that — to avoid inconsistences — one 
must fix an order of integration over the variables £\ and coherently respect it, by defining 
Xi, Zi only in terms of those ty which are integrated before £i in the chosen integration 
order. 

Given a generic cosine diagram (£, the problem of finding for each ti a choice of 
axes Xi, Zi giving rise to an angular measure with the minimal number of residual 
integrations over 9, (f) is — in our knowledge — still unsolved. A partial answer can be 
found in the following proposition, which deals with cosine diagrams of tree type. 

Proposition 1 Given a cosine diagram €. with L Tes vertices and vanishing loop number, 
a choice of axes Xi, Z\ can be given for each ii — in the sense of i f55j) — such that the 
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angular measure is expressed in terms of L rcs residual integrations of cos 8-type ( or, 
equivalently, such that AY has L res elements, all of 9 -type). 

Proof Suppose first that £ is a tree (hence connected by definition). Choose a root 
vertex vi on & Fix arbitrarily a couple of nondegenerate versors (X, Z) and assign 
(xi , Zi ) = (X, Z) for the root vertex and X\ = X for all other vertices v\ of <£. A choice 
of axis Z\ for all vertices vi ^ vi will now be given, by iterating over their distance n 
from the root vertex vi . First define z\ x = £i for all vertices f/ 1 at distance 1 from 
the root — £i being the integration versor associated by definition to the root vertex. 
Then iterate over increasing n > 1 as following: for each vertex v\ n at distance n define 
Zi n+1 = £i n for all vertices vi n+1 departing from vi n and having distance n + 1 from 
the root — £i n being the integration versor associated by definition to Vi n . If £ is not 
connected repeat the above reasoning for each connected component of £, using the fact 
that with this parametrization the angular measure factors in term of angular measures 
over the connected components of the cosine diagram. □ 

In table 0] we will give the complete classification of cosine diagrams for the 
case L re s = 3, as well as the best choice of parametrizations which we were able to 
find according to the lines of reasoning presented in this section (including also some 
additional tricks which will be introduced in the following section 14. 4. 3|) . 

4-4-S- Cosine diagrams with disconnectible lines As mentioned in section 14.4.11 it 
is sometimes useful to encode additional properties of the integrand function F in 
the cosine diagram £, by distinguishing different types of lines. This technique 
comes particularly at hand when parametrizing loop integrations for the amplitudes 
described in equation (|52p. Below we identify a particular situation when the integrand 
function depends on its arguments in a simpler way: an appropriate codification of 
this information in the cosine diagram and a careful analysis may sometimes lead 
to the detection of new types of factored analytically-known angular integrals and, 
consequently, to an additional decrease of the number of residual angular integrations. 
The situation of interest is when the integrand function F can be factored as follows: 

F = (block with no dependence on £ h .£,) ( -= -= — - 

V 1 12 J \l+t 2 h +tl + 2£ h .£ l2 . 



n being a strictly positive integer; in this case the line joining vertices v\ x and Vi 2 in 
£ will be graphically denoted with a dashed style instead of a continuous one, and we 
will say that such a line — together with the associated cosine cos^^ := £i x .£\ 2 between 
loop momenta £^ and £i 2 — is disconnectible. 

The motivation for introducing the notion of disconnectible lines derives from the 
following observation: supposing that we have been able to pick up a parametrization 
of angular measure Yli d(cos9i) d(j>i such that 

(i) cos 9^i 2 belongs to the set of angular integration variables (cos#i, • • • , cos^ res ) 



Automatic parametrization of amplitudes 



43 



(ii) the only dependence of the integrand on cos0^ 2 is via the scalar product i^.i^, 
i.e. no spurious dependencies on cos 9^ of other scalar products are re-introduced 
by the parametrization — see section 14. 4. 2\ 

then the integration 



can be factored away and carried out analytically, thus leading to a lower number of 
required residual integrations. 

However, this is not the actual end of the history: a weaker situation can sometimes 
happen when it is not possible to integrate directly over a disconnectible cosine cos^^, 
but such a cosine can nonetheless be expressed in terms of some of the variables that 
parametrize the angular measure, with the condition that one of these variables must 
appear only in the expression of the scalar product i^.i^. In this case, again, the 
integral over such distinguished angular variable of the power of propagators containing 
£-h-£i 2 — see flUl) — can De factored and, possibly, performed analitically. A typical 
example of such a case is given by the combination 



originating for example in some cases from the "standard" angular parametrization 
defined in (157)1 : it should be noticed that integrals in the form of (|75)) can always be 
performed by means of the technique of complex residues. 

Additional complications arise when we want to minimize the number of angular 
integrations in presence of an integrand function F whose associated cosine diagram (£ 
possesses disconnectible lines. As a matter of fact, in this case we must face two different 
requirements at the same time: as in previous section, we must minimize the number of 
angular integration variables on which F depends; in addition, we must now maximize 
as well the number of occurrences of integrals in the form of (|7H)) and (j7fij1 . 

Such a problem is quite complicated, and finding a solution implies the enumeration 
and the inspection of all possible cosine diagrams with disconnectible lines which can 
arise up to the maximal number of residual loops L Tes we are interested in. 

To give a taste of how things look like, we illustrate in table H] the "simple" case 
L res = 3. In the first column, all possible cosine diagrams are listed; a dashed line means 
that the associated cosine is disconnectible in the sense defined by equation (|74|). As 
explained above, in some cases the corresponding angular measures can be simplified by 
analytically integrating over appropriate angles; thus, in the second column "improved" 
cosine diagrams can be found, describing the cheapest measures which will be used in 
the final parametrization; in this column dashes mean "a 9 has been integrated" — see 
(175)) — , while dots mean "a has been integrated" — see (175)) — . In the third column 
one can see the contribution to the set AV coming from each line of the improved cosine 




(75) 



Jo \l + £ 2 2 + £l + 2e 2 h£2.£3j 

li-£z = cos 6*2 cos 9s + sin #2 sin #3 cos 03, 




(76) 
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Table 4. Our choice of "cheapest" angular measures for the case L ICS — 3. See the 
text for all the explanations about graphical conventions used in this table. 

Initial cosine diagram Final cosine diagram Variables lines depend on Cost 



(0,0) 



(0,0) 



line 1-2 — > {0 2 } 



(1,0) 



(0,0) 



line 1-2 — > {9 2 } 



(1,0) 



line 1-2 
line 1-3 



{0 3 } 



(2,0) 




line 1-2 
line 1-3 



(2,0) 




line 1-2 
line 1-3 



(2,0) 




line 1-2 
line 1-3 



{02} 



(2,0) 





line 1-2 
line 1 3 
line 2 3 



■{02} 
■{03} 

{0 2 ,e 3 ,h} 



(2,1) 



diagram — see (J6U)) . Finally, in the last column we indicate the cost of the final angular 
measure as a couple (number of integrations on 9, number of integrations on (f). 

We just mention that in the more realistic case of L res = 4 we have 65 
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different possible cosine diagrams, which have been analyzed by hand and, successively, 
incorporated in our code. In addition, the program takes automatically care of the fact 
that at L rcs = 4 each basic cosine diagram can manifest in up to 4! = 24 different 
possible permutations. 

5. Implementation, and results 

The techniques described in section 0] were actually used to produce a C++ program; 
given a graphical description of the basis of effective vertices one intends to use — as 
described in section l4~2l — , and a specification of the set of diagrams to be parametrized, 
our code produces 

• an histogram of the distribution of the cost C — as defined in section 14.11 — for 
the specified set of diagrams whose parametrization has been required. Examples 
of such distributions for some different choices of diagrams and bases of effective 
vertices can be found in figure figure HH1 and figure [T7| 

• the C++ functions needed to perform the numerical integration of amplitudes, ready 
to be compiled if an actual numerical implementation of each of the effective 
vertices which have been used during the parametrization stage is given; the set 
of amplitudes is then placed in a dynamical library, with each subroutine ready 
to be retrieved by the integration program. We do not enter here into the very 
complicated details of numerical evaluation of amplitudes, which will be described 
in future publications 

• the specifications needed to compute the RG functions /3 M (<? M ), rj M (g M ), 772m(<?m); 
these informations consist in a list of all required combinatorial and O (N)- factors, 
plus — after a suitable file containing the results of numerical integration has been 
provided — the actual value of each integral in question, with an estimation of its 
numerical error. 

In general, the program presents a lot of useful switches to precisely tune its behaviour: 
for example, subsets of all the needed graphs — or even single graphs — can be selected 
to be parametrized, and a graphical automated output containing all the details of the 
process is provided for debugging purposes if requested; more than one parametrization 
for the same amplitude can be produced and stored when a comparison is desired. 

The code consists of about 14.000 lines of C++ — about 7.000 for graph library and 
other auxiliary libraries, the rest in higher-level routines for the manipulation of scalar 
amplitudes — which have been written and tested from scratch. Being modular in the 
spirit of C++ classes (see [33 ) — with a different class corresponding to each different 
concept exposed in section like substitutions of effective vertices, loop bases, angular 
measures and so on — the code can be quite easily modified to include new computational 
simplifications or new strategies, as discussed in section 14.11 On a typical machine of 
the GHz-era, the 7 loops can be entirely parametrized — and the corresponding C++ code 
produced — in some hours; peak memory occupation is entirely acceptable. 
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1 2 3 4 5 6 

Number of residual integrations 

Figure 15. Total cost of loops from 1 to 5. The histogram accounts for amplitudes 
contributing to v[ 2: °\ dTi°'/dp 2 and rf 4 '°\ including higher-points amplitudes 
obtained from factorization — see section !5.6.2l and section lB.5.11 Used basis of effective 
vertices are 1-loop functions up to the pentagon (see figure 01 plus 2-loop functions of 
figure |3| plus the renormalized sunset. 



Some additional details about the structure, the design and the results of the 
program which we had to omit here due to the lack of space can be found in [4"U] 
(although it documents a slightly outdated version of our framework): for example, the 
explicit parametrization of all amplitudes up to the 4-loop level with "simple" one-loop 
effective vertices can be found there. 

We would now like to discuss in detail some of the typical results produced by the 
program. 

Two histograms detailing the costs of the parametrization of amplitudes for loops 1 
to 5 all together, and for loop 6 alonei:, respectively, are shown in figureUHland figure HE) 
the basis of functions employed here is a "standard" one, composed by the polygons up 
to the pentagon, the renormalized "sunset" effective vertex, all other low-cost 2-loop 
functions shown in figure H3 (see section I3.4.3|) and — from 6 loops and higher — the 
effective vertex corresponding to the 2-loop diagram of figure 0]c. All resulting costs 
are not more than 6; this fact is in agreement with what has been previously found by 
Nickel et al., who declare in j22j that they were able to obtain all their parametrizations 
of 6-loop diagrams with integrations in not more than 6 residual dimensions. 

X By "amplitudes contributing at loop L alone to some collection of quantities X", we mean the set 
of all new amplitudes which have to be evaluated in our I-schcme to obtain all the quantities in X at 
L loop level, assuming that all the amplitudes needed for the evaluation at loop level L — 1 of r' 2 ' '', 
dv[ 2 ' ^ j dp 2 and rf 4 ' '' at vanishing external momenta, are already known. 
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Figure 16. Total cost of the 6-th loop alone. The histogram accounts for amplitudes 



^(2,0) 



and If '°\ including higher-points amplitudes 



contributing to Ii 

obtained from factorization — see section pSE 2l and section l^.5.11 Used basis of effective 
vertices are 1-loop functions up to the pentagon (see figure [HJ plus 2-loop functions of 
figureEl plus diagram of figure0]c, plus the renormalized sunset. 



In spite of the fact that we are working in different renormalization schemes and 
that, consequently, we have a different number of diagrams to evaluate, it would 
be nevertheless quite interesting to be able to perform a more precise quantitative 
comparison between our histograms and the corresponding distributions obtained by 
Nickel and coworkers; however, this information has never been published. Our general 
impression is that our results could be slightly sub-optimal, mainly due to the fact that 
for some "difficult" graphs it is probably possible to work out cheaper parametrizations 
by hand, in the spirit of figure HJ All in all, the results we obtained for the 6-loop case 
can be also interpreted as being very encouraging, if one bears in mind that 

• the simplifications implemented in our code (see the list in section I4.1J1 are just 
a few, and not very complicated; nonetheless, they are enough to obtain a set of 
parametrizations which would seem comparable to that derived in the original work 
of Nickel and coworkers 

• our program is very general, so that other techniques can be easily added if needed 

• a separate treatment of a few particularly expensive graphs is always possible 

• if we want to tackle the computation of the 7-th loop, we must be anyway prepared 
to the appearance of a consistent number of integrals with a residual dimensionality 
much larger than 6. 

Going on to the results relative to the parametrization of the 7-th loop shown 
in figure El — which are the most important since their knowledge represents the 
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0123456789 

Number of residual integrations 




0123456789 

Number of residual integrations 

Figure 17. Total cost of the 7-th loop alone. The histogram accounts for amplitudes 
contributing to v[ 2 '°\ dTi'°^/dp 2 and r| 4 '°\ including higher-points amplitudes 
obtained from factorization — see section 13.6.21 and section 13.5.11 Used bases of 
effective vertices are: above, 1-loop functions up to the pentagon (see figure|BJ) plus the 
renormalized sunset; below, 1-loop functions up to the pentagon, plus 2-loop functions 
of figure |31 plus diagram of figure 01 c, plus the renormalized sunset. 

motivation for this work, in view of their possible use to improve available estimates 
of critical quantities — the situation looks even more promising. The larger residual 
dimensionality is here 9, still low enough to make possible a high-precision numerical 
evaluation of the amplitudes if enough computer power is available. By comparing the 
two profiles shown in the figure one can notice how the two-loop effective functions 
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play an essential role in lowering the overall complexity of the computation of a 
factor of ~ 3 (this rough estimate can immediately be obtained by remembering, as 
explained in section 13 .3[ that the main part of computer time will be spent in the 
evaluation of the integrals of highest dimensionality). Such results bring us to the 
remarkable conclusion that even at a 7-loop level the approach of Baker, Nickel, Green 
and Meiron remains fruitful, leading — within the present framework — to an acceptable 
set of parametrizations for all the required diagrams, including those contributing to the 
4-point function. 

As a final remark, some preliminary looks at the corresponding results for the 2- 
point functions at 8 loops (not reported in this article) seem to indicate that the set of 
effective vertices which allows a successful parametrization of the 7 loops would not be 
adequate to conveniently parametrize the 8 loops; anyway, due to the number and the 
complexity of required integrals the possibility of carrying out the computation of the 
8-th loop is purely academic at the moment. 

6. Epilogue 

In this work we present a solution to the problem of building an automated framework for 
the parametrization of Feynman amplitudes needed to compute — at the level of 7 loops 
and higher — the Renormalization-Group /3-function and the anomalous dimensions r], r\2 
for the 3-dimensional O(N) vector model. 

The core of the employed computational method — in the spirit of [1 -|2J — consists 
in replacing effective-vertex functions and applying other analytic simplifications to the 
original expressions for amplitudes in momentum space. 

Some of the key features of our own approach, as compared to that in PQ-|2], are 
the following: 

• a different intermediate renormalization scheme — the I-scheme — has been chosen, 
which has simpler renormalization properties as well as a Callan-Symanzik equation 
with exactly-known coefficients, see section l2~2l 

• a systematic treatement of the parametrization of loop integrations in momentum 
space (leading in our framework to a "good choice of angular measures" ) has been 
developed in section |U 

• an extended set of effective vertices, including some convenient 2-loop functions, 
has been used to perform the analysis. 

The employed technique can eventually allow to reformulate the needed amplitudes in 
terms of relatively low-dimensional integral representations; however, a general proof 
of the effectiveness of the method is missing, and its practical applicability to the full 
7-loop case had not been demonstrated so far. 

In spite of the theoretical possibility that some topological obstructions could show 
up, leading to some very difficult graph with a very high number of residual integrations, 
the automated parametrization of all amplitudes needed to compute RG-functions 
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revealed that — luckily — at the level of 7 loops nothing similar happens: this work shows 
that a reasonably small set of computational simplifications is still sufficient to obtain 
for all the needed amplitudes efficient parametrizations with a residual dimensionality 
of not more than 9, thus opening the way to a full 7-loop evaluation of RG /^-function. 

On the other hand, the large number of needed diagrams (~ 4000) and the necessity 
of carrying out each integration with 7-8 digits of precision in dimensionalities up to 9 
nonetheless classify the problem as a formidable one in terms of required computational 
resources. How to solve it is an open issue, and possible solutions will be discussed in 
an incoming publication. 
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Appendix A. Renormalization of the O(N) vector model and 
Callan-Symanzik equations 

After a brief review of our notations, we give in this section a concise presentation of 
the O(N) vector model, of its renormalization and of Callan-Symanzik equations. The 
reader looking for a more detailed and pedagogical introduction should refer e.g. to [T8] . 

a m, nzi- 

We will work in the Euclidean space M 3 , and consider an 0(AQ-invariant quantum 
field theory built out of scalar fields <fi l {x) in a fundamental representation of O(N) - 
where i = 1 . . . , N. 

For shortness we will omit field arguments, O(N) or M. 3 indices unless they play a 
crucial role; we will write x instead of a collection of objects (xj G §>)™ =1 if no ambiguity 
arises (that is, if the set § and the value n can be clearly deduced from the context). In 
this notation, the Fourier transform of a function f(xi, . . . , x n ) invariant under global 
space translations is written 

/(g) (2vr) 3 6 3 (J2r) ■= J dx exp(ip.x) f(x) (A.l) 

meaning as usual 

n „ / n \ 

f( Pl , . . . , pj (2vr) 3 5 3 QT>,) == J f IP 3 ^ exp^.a^J f(x u . . . , x n ). (A.2) 
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The bare action with sources J B , K B is defined by: 

S B (0 B , Jb, K b ) := J d 3 x Q(V0 B ) 2 + ^0 B + ^(0 B ) 2 - J B B + \K»<£) . (A.3) 
It should be noticed that 

(^*b) = ^hhhu^B^B^B^B ^hhkiU - = g (^1112^3*4 ^ii 13^214 ^hU^hia) • (A-4) 

The generating functional for the bare theory is defined by: 

Z B (J B ,K B ):= y'[rf0 B ]exp(- 1 S B (0 B , J B , K B )) (A.5) 

(an appropriate regularization parametrized by an ultraviolet cutoff A is assumed, to 
give meaning to functional integration). The effective action T B is defined via Legendre 
transform of the generator of connected diagrams, logZ B : 

logZ B (J B ,K B ) +r B (0 B ,^ B ) := J d 3 x J B (x)0 B (x) b := ^j^- (A.6) 

One-particle irreducible (OPI) correlators are defined as usual from the expansion of the 
effective action w.r.t. its sources: 



°° I f 

b\v>b,K b ) = ~fiw~\ j d xdyr B E ' V2 \x-,y) B (xi) . . . <p B (x E ) tf B (l/i) • • • K B (y V2 ). 

E,V 2 =0 ' 2 ' J 

Exploiting 0(A r )-invariance of the theory, we factor some correlators evaluated at zero 
momenta into a 0(AQ-invariant tensor times a 0(A^)-scalar function: 

r B 2 ' 0) n i2 (p 1 ,P 2 ; • • 0Ie=o == S ili2 lf' 0) (- • •) (A.7) 
rf'\ , 2 (Pi , p 2 ; g; ■ ■ •) | E=9 =o =: S h n r B 2 ' 1} (• • •) (A.8) 

Tb ' ^iii 2 i3i 4 (Pl)P2)P3)P4i ' ' ')lp=0 =: ^ii i 2 13 M ' ''(''') (A-9) 

(the dots • • • in expressions above denote the dependence on other parameters) . 

In a generic renormalization scheme, bare quantities are related to corresponding 
renormalized ones according to the following transformations: 



B =:v^0 R (A. 10) 

m 2 B =: (ml + 8ml) / (A.ll) 

A B =: A R (A.12) 

^[0 2 ] r ,:=^^0b + <)- (A. 13) 

The transformation from bare to renormalized quantities is performed at fixed 
regularization parameter A and is parametrized by the scheme-dependent quantities 
Sml and Z's. The dependence on the renormalization scheme will be emphasized 
by capitalized suffixes such as m or i. The special suffix r is used to label scheme- 
dependent quantities in relations holding for a generic renormalization scheme with non- 
zero renormalized mass. The same conventions used in the case of bare quantities and 
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their Legendre and Fourier transforms will be adopted also for renormalized quantities 
(apart from a straightforward change of suffix). 

The renormalized (sourceless) action is obtained from the bare one after expressing 
it in terms of renormalized quantities: 



SM ■= S B [4> B = y/Z^^Jn = 0,K B = Oj (A.14) 

= j d 3 x (%^(V0 R ) 2 + \ (ml + 5m R ) R + ^(0 R ) 2 ) . (A.15) 

The renormalized functional generator Z R , the renormalized effective action T R and OPI 
correlators are expressed in terms of bare quantities as follows: 

Z R {J R ,K R ) := J[d<t> n ] exp(-S R (0 R ) + J J R R - J ^ R [0 2 ] R J (A.16) 



,2 



Z*\-h =^=,K B = ^K R )e W (~ j l\"r ; Z\ Ti 

2 \ 

1 



r R (0 R ,ir R ) = T B (0 B = y/Z^Kv = -^-fK n j +-j K R Zf; R (A. 17) 

rf^(*; y) = {Z^ (^A Y^\x; y) + \ 5 E0 5y 2 1 <. (A.18) 

In the rest of this article we will only deal with correlators with E > 0, so we do 
not consider anymore the local contact term in the right hand side of (|A.18|) : however, 
it should be noticed that in the case (E, V 2 ) = (0, 1) the contact term does modify 
Callan-Symanzik equations (|A.21jl . 

The Callan-Symanzik differential operator in a generic scheme is defined as 

(A. 10) 



V R := m R - 



R A B ,A 



Callan-Symanzik equations are obtained by application of the Callan-Symanzik operator 
on both sides of (Fourier transform) of (jA.18|) and by use of identities 

(p;q) = V K [ml\ feq) = P R [m B ] Y? > W) ( B q,q V2+1 = 0) (A.20) 

which follow from the chain rule and from the fact that a derivative of bare OPI 
correlators w.r.t. the bare mass is equivalent — in non-pathological cases — to an 
insertion of B /2 at zero momentum. 

The resulting Callan-Symanzik equations in a generic renormalization scheme are 
(assuming E > 0): 

d „ , , d E 



m 



T?™ (p;g) 



dm R dg R 2 

= ml a R (g R ) if > V2+1) (p; q, q V2+1 = 0) (A.21) 
where we introduced the adimensional renormalized coupling 

g R := A R /m R (A.22) 
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and the RG-functions 



V R [g R ] 
V R [\ogZ^] 



(A.23) 
(A.24) 
(A.25) 



0"r(^r) : 




(A.26) 
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